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Abstract 


Industrial  assembly  involves  sensing  the  pose  (orientation  and  position)  of  a  part. 
Efficient  and  reliable  sensing  strategies  can  be  developed  for  an  assembly  task  if  the 
shape  of  the  part  is  known  in  advance.  In  this  paper  we  investigate  two  problems  of 
determining  the  pose  of  a  polygonal  part  of  known  shape.  In  the  first  problem,  the 
part  has  a  continuum  of  possible  poses,  while  in  the  second  problem,  it  has  a  finite 
number  of  possible  poses. 

More  specifically,  the  first  problem  involves  determining  the  pose  of  a  convex  n-gon 
from  a  set  of  m  supporting  cones,  that  is,  cones  with  both  sides  supporting  the  polygon. 
An  algorithm  with  running  time  0(nm)  which  almost  always  reduces  to  0(n  +  m  log  n) 
is  presented  to  solve  for  all  possible  poses  of  the  polygon.  As  a  consequence,  the  polygon 
inscription  problem  of  finding  all  possible  poses  for  a  convex  n-gon  inscribed  in  another 
convex  m-gon,  can  be  solved  within  the  same  eisymptotic  time  bound.  We  prove  that 
the  number  of  possible  poses  cannot  exceed  6n,  given  m  >  2  supporting  cones  with 
distinct  vertices.  Experiments  demonstrate  that  two  supporting  cones  are  sufficient  to 
determine  the  real  pose  of  the  n-gon  in  most  cases.  Our  results  imply  that  sensing  in 
practice  can  be  carried  out  by  obtaining  viewing  angles  of  a  planar  part  at  multiple 
exterior  sites  in  the  plane.  As  a  conclusion,  we  generalize  this  and  other  sensing 
methods  into  a  scheme  named  sensing  by  inscription. 

On  many  occasions  a  parts  feeder  will  have  reduced  the  number  of  possible  poses 
of  a  part  to  a  small  finite  set.  In  order  to  distinguish  between  the  remaining  poses  of 
the  part  some  simple  sensing  or  probing  operation  may  be  used.  Our  second  problem, 
called  sensing  by  point  sampling,  is  concerned  with  finding  the  minimum  number  of 
sensing  points  required  to  distinguish  between  a  finite  set  of  polygonal  shapes.  In 
practice  this  can  be  carried  out  by  embedding  a  series  of  point  light  detectors  in  a 
feeder  tray  or  by  using  a  set  of  mechanical  probes  that  touch  the  feeder  at  a  finite 
number  of  predetermined  points. 

Intuitively,  each  sensing  point  can  be  regarded  as  a  binary  bit  that  has  two  values 
‘contained’  and  ‘not  contained’.  So  the  robot  senses  a  shape  by  reading  out  the  binary 
representation  of  the  shape,  that  is,  by  checking  which  points  are  contained  in  the 
shape  and  which  are  not.  The  formalized  sensing  problem:  Given  n  polygons  with  a 
total  of  m  edges  in  the  plane,  locate  the  fewest  points  such  that  each  polygon  contains 
a  distinct  subset  of  points  in  its  interior.  We  show  that  this  problem  is  equivalent  to 
an  NP-complete  set-theoretic  problem  introduced  as  Discriminating  Set.  By  a  reduc¬ 
tion  to  Hitting  Set  (and  hence  to  Set  Covering),  an  0{n’^m'^)  approximation  algorithm 
is  presented  to  solve  the  sensing  problem  with  a  ratio  of  2 Inn.  Based  on  a  reverse 
reduction,  we  prove  that  one  can  use  an  algorithm  for  Discriminating  Set  with  ratio 
clogn  to  construct  an  algorithm  for  Set  Covering  with  ratio  clogn-l-O(loglogn).  Thus 
approximating  Discriminating  Set  exhibits  the  same  hardness  as  that  of  approximating 
Set  Covering  recently  shown  in  [44]  and  [7];  this  result  implies  that  the  ratio  2 Inn  is 
asymptotically  optimal  unless  NP  C  DTIME(nP‘’*’'*°*’*).  We  also  analyze  the  complex¬ 
ity  of  subproblems  of  Discriminating  Set,  based  on  their  relationship  to  a  generalization 
of  Independent  Set  called  3-Independent  Set.  Finally,  we  present  experimental  results 
showing  that  sensing  by  point  sampling  is  mostly  applicable  when  poses  are  densely 
distributed  in  the  plane. 


(D  (b) 

Figure  1:  Geometric  sensing  of  an  8-gon.  (a)  The  case  of  a  continuum  of  possible  poses:  Four 
incidence  constraints  and  I2  suffice  to  determine  the  pose,  (b)  The  case  of  a  finite  number 

(three)  of  possible  poses:  Only  the  containment  of  two  points  qi  and  92  needs  to  be  checked  to 
determine  the  real  pose. 

1  Introduction 

Sensing  a  part  involves  determining  both  its  shape  and  pose.  By  pose  we  mean  the  position 
as  well  as  the  orientation  of  the  part.  Prior  to  selecting  a  sensing  method,  we  often  will  make 
some  assumptions  about  the  shape  of  the  part  to  be  sensed.  The  resulting  sensing  method 
is  affected  greatly  by  what  is  known  about  the  shape.  For  instance,  without  making  any 
assumptions,  we  might  not  even  be  able  to  start  segmentation  of  the  part  image,  whereas 
knowing  that  the  shape  is  convex  polygonal,  we  can  employ  some  simple  non- vision  technique 
such  as  finger  probing.  An  effective  sensing  method  should  m2dce  use  of  its  knowledge  about 
the  part  shape  os  much  as  possible  to  attain  simplicity,  efficiency  and  robustness. 

Parts  in  many  assembly  applications  are  manufactured  to  high  precisions,  so  we  can  make 
the  assumption  that  their  shapes  are  known  reasonably  well  in  advance.  Accordingly,  the 
design  of  sensing  strategies  should  be  based  on  the  geometry  of  parts.  The  task  of  sensing 
reduces  to  obtaining  enough  geometric  constraints  which,  when  combined  with  the  part 
geometry,  suffice  to  derive  the  part  poses.  Consequently  minimizing  the  necessary  geometric 
constraints  becomes  very  important  for  reducing  the  sensing  complexity.  In  this  paper,  we 
propose  two  approaches  for  sensing  polygonal  parts  in  known  shapes,  one  applicable  to  a 
continuum  of  possible  poses  and  the  other  applicable  to  finite  possible  poses. 

Perhaps  the  simplest  geometric  constraint  on  a  polygon  is  incidence — when  some  edge 
touches  a  fixed  point  or  some  vertex  is  on  a  fixed  line.  For  instance.  Figure  1(a)  shows 
an  8-gon  constrained  by  two  points  pi,p2  and  two  lines  /i,/2.  The  question  we  want  to 
ask  is:  Generally  how  many  such  constraints  are  necessary  to  fix  the  polygon  in  its  real 
position.  Note  that  any  two  such  incidence  constraints  will  confine  all  possible  positions  to  a 
locus  curve  which  consists  of  a  finite  number  of  algebraic  curves  parameterized  by  the  part’s 
orientation.  Three  constraints,  as  long  as  not  defined  by  collinear  points  or  concurrent  lines, 
will  allow  only  a  finite  number  of  valid  poses.  These  poses  occur  when  different  locus  curves, 
each  given  by  a  pair  of  constraints,  intersect  at  the  same  orientations.  An  upper  bound  on 
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the  number  of  possible  poses  can  be  analyzed.  The  addition  of  a  fourth  constraint  is  usually 
enough  to  reduce  this  number  to  one — the  real  pose. 

If  all  incidences  are  given  by  lines,  sensing  can  be  viewed  as  inscribing  the  polygon  into 
a  larger  polygon  (not  necessarily  bounded)  defined  by  these  lines;  if  all  constraints  are  given 
by  points,  sensing  can  be  viewed  as  placing  the  polygon  defined  by  these  points  into  the 
sensed  polygon. 

Point  constraints  can  be  created  by  “probing”  the  polygon  along  various  directions  with 
a  tactile  sensor  or  a  range  finder.  Line  constraints  can  be  obtained  with  an  angular  sensor 
scanning  across  the  object  at  exterior  sites.  In  Section  2  we  study  the  case  of  sensing  with 
line  constraints  only,  offering  a  very  efficient  algorithm  to  solve  for  all  possible  poses.  We 
give  a  tight  upper  bound  on  the  number  of  poses  for  three  lines,  and  conduct  experiments 
to  show  that  four  lines  are  practically  enough  to  determine  the  real  pose. 

The  set  of  possible  poses  can  often  be  reduced  from  a  continuum  to  a  finite  number  in 
advance,  by  planned  manipulations  such  as  pushing  or  squeeze-grasping,  or  by  sensing  from 
geometric  constraints  as  mentioned  above.  More  specialized  sensing  methods  can  be  devised 
to  distinguish  between  the  remaining  finite  number  of  poses.  We  now  view  each  pose  as  a 
closed  set  of  points  in  the  plane  occupied  by  the  part  in  that  pose,  so  that  sensing  becomes 
distinguishing  point  sets,  say  . . .  ,P„.  An  easy  way  is  to  sample  several  points,  checking 
which  of  them  are  contained  by  the  current  point  set  (pose).  Suppose  the  same  8-gon  is 
known  to  have  only  three  possible  poses,  as  shown  in  Figure  1(b),  then  the  real  pose  (drawn 
in  solid  line)  can  be  determined  by  verifying  that  it  contains  both  points  qi  and  92-  This 
can  be  implemented  in  a  number  of  different  ways,  for  instance,  by  placing  light  detectors 
underneath  the  point  locations  or  by  probing  the  point  locations  from  above  with  a  robot 
finger. 

In  Section  3  we  address  the  problem  of  distinguishing  finite  polygons  by  point  sampling. 
We  prove  that  minimizing  the  total  number  of  sampling  points  is  NP-complete,  and  offer 
an  approximation  algorithm  with  a  greedy  heuristics.  The  algorithm  produces  a  set  of 
stimpling  points  whose  size  is  within  a  factor  of  21ogn  of  the  optimal.  We  also  exhibit  a 
proof  demonstrating  the  hardness  of  improving  this  approximation  ratio. 


2  Sensing  by  Inscription 

2.1  The  Inscription  Problem 

In  this  section  we  will  study  the  problem  of  detecting  the  pose  of  a  convex  polygon  in  the 
plane  by  taking  views  of  the  polygon  from  multiple  exterior  sites.  The  shape  of  the  polygon 
is  assumed  to  be  known  in  advance,  but  the  pose  of  the  polygon  can  be  arbitrary.  Each  view 
results  in  a  cone  formed  by  the  two  outermost  occluding  rays  starting  from  the  viewing  site. 
The  cone  in  turn  imposes  a  constraint  on  the  possible  poses  of  the  polygon — the  polygon 
must  be  contained  in  the  cone  and  make  contact  with  both  its  sides.  A  containment  in  which 
every  edge  of  the  containing  object  touches  the  contained  object  is  called  an  inscription,  so  we 
shall  say  that  the  polygon  is  inscribed  in  the  cone.  Such  constraints  imposed  by  individual 
views  together  allow  only  a  small  number  of  possible  poses  of  the  polygon,  which  often 
reduces  to  one.  For  example,  Figure  2  illustrates  two  views  taken  of  a  convex  6-gon  P  in 
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Figure  2:  Sensing  the  pose  of  a  polygon  by  taking  two  views. 


some  unknown  pose  from  sites  Oi  and  02  respectively:  The  two  cones  Ci  and  C2  thus  formed 
determine  the  real  pose  of  P,  and  this  pose  can  be  solved  using  the  algorithm  presented  later 
in  Section  2.2.4. 

The  above  sensing  approach  appears  to  be  simple,  but  to  make  it  efficient  and  to  minimize 
the  cost  of  sensing  hardware  we  would  like  to  take  as  few  views  as  possible.  This  leads  us  to 
the  main  question  of  this  section:  How  many  views  are  sufficient  in  the  general  case  in  order 
to  determine  the  pose  of  a  convex  n-gon?^ 

The  answer  to  the  above  question  is  two,  and  to  argue  this  answer  we  will  go  through 
several  steps  each  of  which  occupies  a  separate  subsection:  Section  2.2  describes  how  to 
compute  the  set  of  possible  poses  for  a  convex  polygon  inscribed  in  multiple  cones,  and 
derives  an  upper  bound  on  the  number  of  possible  poses  for  two-cone  inscription  in  particular; 
Section  2.3  empirically  demonstrates  that,  in  spite  of  the  upper  bound,  two  cones  turn  out  to 
be  suflBcient  in  most  cases  to  uniquely  determine  the  pose  of  an  inscribed  polygon.  Section  4 
further  discusses  the  extensions  of  this  method  and  proposes  a  general  sensing  scheme — 
sensing  by  inscription. 

2.1.1  Related  Work 

Canny  and  Goldberg  [11]  have  introduced  a  Reduced  Intricacy  in  Sensing  and  Control  (RISC) 
paradigm  that  aims  at  improving  the  accuracy,  speed  and  robustness  of  sensing  by  coupling 
simple  and  specialized  hardware  with  fast  and  task-oriented  geometric  algorithms. 

The  cross-beam  sensing  method  developed  in  [55]  finds  the  orientation  of  a  polygon 
(or  polyhedron)  by  measuring  its  diameters  along  three  different  directions  and  comparing 
the  measurements  with  the  precomputed  diameter  function  [27];  then  it  solves  a  vertex¬ 
line  correspondence  problem  for  the  position  of  the  polygon  by  least  squares  fitting.  This 

should  be  noted  that  there  exist  cases  in  which  the  pose  of  a  convex  n-gon  cannot  be  uniquely 
determined  no  matter  how  many  views  are  taken.  This  happens  only  if  the  polygon  preserves  self-congruence 
over  certain  rotations.  (It  is  not  hard  to  see  that  in  such  a  case  the  rotation  angle  must  be  a  multiple  of 

where  ^  is  a  positive  integer  such  that  ib|n.)  However,  all  congruent  poses  are  usually  considered  as  the 
same  in  the  real  applications. 


method  essentially  determines  the  pose  by  inscribing  the  polygon  in  a  hexagon  constructed 
from  the  sensory  data.  The  idea  of  characterizing  shapes  with  diameters  and  chords  was 
also  addressed  earlier  in  [53]. 

For  the  special  case  that  the  poses  are  finite,  [49]  presents  an  efficient  method  of  placing 
a  registration  mark  on  the  object  so  that  the  pose  can  be  recognized  by  locating  the  mark 
position  (with  a  simple  vision  system).  For  robustness  to  sensor  imperfections,  the  marked 
point  maximizes  the  distance  between  the  nearest  pair  among  its  possible  locations.  For  the 
case  that  the  number  of  parts  is  finite,  [28]  recognizes  a  part  with  a  modified  parallel-jaw 
gripper  by  a  sequence  of  grasp  actions  and  diameter  measurements. 

Model-baised  recognition  and  localization  can  often  be  regarded  as  a  constraint  satis¬ 
faction  problem  which  searches  for  a  consistent  matching  between  sensory  data  (e.g.,  2-D) 
and  model(s)  (e.g.,  3-D)  based  on  the  geometric  constraints  between  them.  (See  [30].)  [26] 
discusses  how  to  identify  and  locate  a  polyhedron  on  a  known  plane  using  local  information 
from  tactile  sensors  which  includes  the  position  of  contact  points  and  the  ranges  of  surface 
normals  at  these  points.  Motivated  by  an  interpretation  tree  developed  in  [26],  [52]  deter¬ 
mines  the  pose  of  an  object  grasped  by  a  hand,  under  a  situation  very  close  to  inscription. 
Also  using  an  interpretation  tree  search,  [37]  solves  the  pose  of  a  polyhedra  by  matching  a 
set  of  3-D  line  segments,  obtained  by  three  light-stripe  range  finders,  to  model  faces;  then 
the  pose  uncertainty  is  estimated  using  the  covariance  matrix  of  the  endpoints  of  these  line 
segments.  [14]  uses  a  polynomial  approach  to  solve  for  the  line-to-plane  correspondences 
involved  in  pose  determination.  Based  on  a  generalized  Hough  transform,  [42]  estimates  the 
pose  of  a  3-D  object  by  matching  point  triples  on  the  object  to  possibly  corresponding  triples 
in  the  image. 

In  the  meantime,  a  variety  of  polygon  shape  descriptors  [2],  [45]  have  been  analyzed 
theoretically  and/or  demonstrated  experimentally  to  be  efficient  and  robust  to  uncertainties. 

Geometric  algorithms  for  sensing  unknown  poses  as  well  as  unknown  shapes  have  also 
been  studied.  Cole  and  Yap  [16]  considered  “finger”  probing  a  convex  n-gon  (n  unknown) 
along  directed  lines,  and  gave  a  procedure  guaranteed  to  determine  the  n-gon  with  3n  probes. 
This  work  Wcis  later  extended  in  [18]  which  investigates  the  complexity  of  various  models 
for  probing  convex  polytopes  in  d-dimensional  Euclidean  space.  Using  a  more  powerful 
probe  model  that  returns  not  only  the  contact  points  but  also  the  normals  at  these  points, 
Boissonnat  and  Yvinec  [9]  showed  how  to  compute  the  exact  shape  of  a  simple  polygon  as 
long  as  some  mild  conditions  about  the  shape  are  met.  [40]  gives  algorithms  of  projection 
probing  and  line  probing  (similar  to  inscription)  which  perform  the  numbers  of  probes  that 
are  both  sufficient  and  necessary  to  determine  a  convex  n-gon.  Close  upper  and  lower  bounds 
are  derived  in  [41]  on  the  number  of  composite  probes  to  reconstruct  a  convex  n-gon,  where 
a  composite  probe  comprises  in  parallel  several  supporting  line  probes  or  finger  probes. 

The  polygon  containment  problem,  that  is,  deciding  whether  an  n-gon  P  can  fit  into  an 
m-gon  Q  under  translations  and/or  rotations,  has  been  studied  by  various  researchers  in 
computational  geometry.  (See  [5],  [6],  [12],  [23].)  In  the  case  where  Q  is  convex,  the  best 
known  algorithm  runs  in  time  O(m^n)  when  both  translations  and  rotations  axe  allowed  [12]. 
Here  we  will  deal  with  a  special  case  of  containment  in  which  each  edge  of  Q  must  touch 
P]  this  constraint  causes  a  reduction  of  the  running  time  to  0{mn),  or  0{n  -1-  mlogn)  in 
practical  situations. 
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(a)  (b) 

Figure  3:  A  triangle  sliding  in  a  cone.  Vertices  pi  and  p2  move  along  the  two  sides  of  cone  C. 
The  locus  of  po  is  an  elliptic  curve  (possibly  degenerating  into  a  line  segment)  parameterized  by  the 
angle  7  between  directed  edge  pjpi  and  the  x  axis.  There  are  two  different  c^ses:  (a)  po  is  above 
edge  P1P2;  (b)  po  is  below  edge  pip2. 

2.2  Multi-Cone  Inscription 

To  simplify  the  presentation,  let  us  agree  throughout  Section  2.2  that  all  angles  take  values 
in  the  half-open  interval  [0, 27r).  In  accordance  with  this  agreement,  any  expression  E  on 
angles  equated  with  or  assigned  to  an  angle  0  in  formulas  such  dsO  —  E  or  9*—  E  will 
be  regarded  as  '‘E  mod  27r’,  even  though  we  do  not  mention  so  explicitly.  (However,  this 
does  not  apply  to  conditions  such  as  E\  =  E2  where  both  sides  are  expressions  on  angles.) 
Moreover,  intervals  for  .angle  values  whose  left  endpoints  are  greater  than  right  endpoints 
are  allowed;  for  example,  an  interval  [q,0],  where  0  <  <  o  <  2x,  is  understood  as  the 

interval  union  [a,27r)  U  [0,y9]. 


2.2.1  A  Triangle  Sliding  in  a  Cone 

We  first  deal  with  the  case  of  a  triangle  in  a  cone,  not  only  because  it  is  the  simplest  but 
also  because  the  case  of  a  polygon,  as  we  will  see  later,  can  be  decomposed  into  subcases 
of  triangles.  Let  ApoPiPz  be  a  triangle  inscribed  in  an  upright  cone  C  with  angle  ^  and 
vertex  o,  where  0  <  <f>  <  ir,  making  contacts  with  both  sides  of  the  cone  at  vertices  pi  and 
P2  respectively.  What  is  the  locus  of  vertex  po  as  edge  pip2  slides  against  the  two  sides  of 
the  cone? 

Two  different  situations  can  occur  with  this  inscription:  (a)  po  is  outside  Apiopa,  and  (b) 
Po  is  inside  (only  when  ZpipoP2  >  <f>)-  (See  Figure  3.)  Assume  that  ApopiP2  may  degenerate 
into  any  one  of  its  edges  but  not  a  point;  writing  a  =  ZpoP2Pi,  ^  =  ZpoPiP2,  r  =  IpoPzI  and 
5  =  \P\P2\1  this  degeneracy  is  taken  into  account  by  the  following  constraints: 


0<a,/9<7r,  s  >  0,  and 


{r  >  0,  if  s  =  0  or  a  >  0; 
0  <  r  <  s,  otherwise. 


Let  us  set  up  a  coordinate  system  with  the  origin  at  o  and  the  y  axis  bisecting  angle  as 
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shown  in  Figure  3.  Then  the  orientation  of  ApoPiP2  can  be  denoted  by  the  angle  7  between 
vector  p2pi  and  the  x  axis.  The  range  of  valid  7  values  can  be  easilj  determined  as 


.V  6  T  <f>  .  .3x  (f>  Zt:  6 


in  case  (a),  and  as  '^  =  'f'  mod  27r  in  case  (b),  where 


3ir  S 

T  +  l  + 


<  27r  +  -  -  -  - 
“  2  2^ 


For  any  valid  7,  there  exists  a  unique  pose  of  the  triangle  in  cone  C;  this  allows  us  to 
parameterize  the  locus  (x,  y)  of  po  hy  7. 

The  two  cases  (a)  and  (b)  yield  results  that  resemble  too  much  to  each  other  to  be  treated 
separately  here.  We  begin  with  writing  out  the  following  equations  for  the  locus  of  po: 


X  =  r  cos(7  —  a)  ±  |op2 1  sin 

£1 

y  =  rsin(7-a)  + |op2|cos^, 


where  notation  ‘±’  means  ‘+’  in  case  (a)  and  ’  in  case  (b) 

and  the  notation  means  just  the  opposite.  Severed  steps  of  manipulation  on  the  above 
equations  plus  a  detailed  subsequent  analysis  on  the  ranges  of  angles  reveal  the  locus  of  po 
as  described  below. 

Namely,  as  edge  piP2  slides  in  the  cone,  po  traces  out  an  elliptic  curve  C  with  implicit 
equation 

ax^  ±  hxy  +  cy^  =  d, 

where 


a  = 


6  = 


2  sin(fq:o)  , 
r  —  rs — - h 


sm 


I  2-2cos^’ 


(2r  cos  a  —  s)s 
sin^  ’ 


r  —  rs 


cos(f  To) 


+  :r- 


cos 


I  2  +  2cos^’ 


d  = 


sin(^T^x\2 


Furthermore,  if  the  orientation  7  changes  monotonically  within  its  range,  po  moves  mono- 
tonically  along  C  except  when  C  degenerates  into  a  line  segment.  In  that  degenerate  case,  po 
moves  along  a  segment  of  a  line  through  the  the  cone  vertex  o  and  with  equation 


{cos  I  X  q:  sin  |  y  =  0, 

cos(|  ^  a)  X  ±  sin(|  T  «)  y  =  0, 


if  r  =  0  (and  thus  a  =  0); 
if  r  ,4  0,  s  0,  and  ;  = 
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Figure  4:  A  convex  polygon  P  rotating  in  a  cone.  The  pose  of  P  is  denoted  by  the  position  of  vertex 
Po  and  the  angle  0  between  directed  edge  poPi  a-nd  the  x  axis.  The  space  of  orientations  [0,2]r)  is 
partitioned  into  closed  intervals,  each  defining  an  elliptic  curve  that  describes  the  corresponding 
locus  of  pq. 

crossing  the  same  point  at  any  two  valid  orientations  71  ^  72  with  71  +  72  =  2ir  ±  ^  when 

the  line  segment  assumes  the  first  equation,  or  with  71  +  72  =  27r  —  ^  +  2a  in  case  (a)  and 

with  7i  +  72  =  2a  +  ^  in  case  (b)  when  it  assumes  the  second  equation. 

Interestingly,  we  observe  that  in  case  (a)  if  r  0  and  j  ,  the  second  line 

equation  reduces  to  a:  =  0  when  ^  =  2a  in  case  (a);  but  this  reduction  does  not  occur  in  case 

(b)  because  j  implies  that  0  <  ^  +  2a  <  27r,  thereby  establishing  sin(|  +  a)  >  0. 

This  observation  reflects  a  small  asymmetry  between  the  two  cases. 

1 

2.2.2  One- Cone  Inscription 

Now  consider  the  case  that  a  convex  n-gon  P  with  vertices  po,pi, . . .  ,Pn-i  in  counterclockwise 
order  is  inscribed  in  a  cone  C.  Let  us  choose  the  seime  coordinate  system  as  in  Section  2.2.1. 
Then  the  pose  of  P  is  uniquely  denoted  by  the  locus  of  some  vertex,  say  po,  and  the  angle  6 
between  the  x  axis  and  some  directed  edge,  say  poPi-  Cleaurly  any  orientation  6  gives  rise  to  a 
unique  pose  of  P;  so  we  can  compute  the  locus  of  po  as  a  function  of  0  over  [0, 27r).  Let  pi  and 
Pr  be  the  two  vertices  currently  incident  on  the  left  and  right  sides  of  cone  C  respectively. 
(See  Figure  4.)  As  long  as  p/  eind  pr  remain  incident  on  these  two  sides  respectively,  the 
problem  reduces  to  the  case  of  Apop/Pr  sliding  in  cone  C  except  that  the  locus  of  po  (an 
elliptic  curve)  now  needs  to  be  parameterized  by  instead  of  7  which  we  had  used  before. 
This  is  easy  for  we  observe  that,  when  p/  ^  p,.  ^  po, 

{9  +  ^PxPqPt  +  ^PoPtPi  —  If,  if  Po  is  above  edge  pipr] 

^  +  ^PiPoPi  +  ^PoPiPr  —  ’T,  otherwise. 

The  three  special  cases  p/  =  po,  Pr  =  Po  and  pi  =  pr  can  be  hemdled  by  substituting  tt  for 
ZpiPoPr  and  Zpipopr  for  ZpoPrP/  for  the  first  case,  tt  for  LpiPopi  and  Ipipopi  for  LpopiPr  for 
the  second  case,  and  0  for  Lpoprpi  (or  IpoPiPr)  for  the  third  case. 
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More  observations  show  that  the  entire  range  [0,27r)  of  orientations  can  be  partitioned 
into  a  sequence  of  closed  intervals,  within  each  of  which  the  vertices  p/  and  pr  incident  on 
cone  C  are  inv2uriant. 

We  present  a  linear-time  algorithm  that  computes  the  above  orientation  intervals  as  well 
as  the  corresponding  elliptic  curves  describing  the  locus  of  po-  The  algorithm  rotates  the 
polygon  counterclockwise  in  the  cone,  generating  a  new  interval  whenever  one  (or  both)  of 
the  incident  vertices  p/  and  p,  changes;  the  new  incident  vertex  (or  vertices)  is  determined  by 
a  comparison  between  angle  ^  and  the  angle  intersected  by  the  two  rays  pi-ipi  and  prPr-i- 
Let  [^niinj  ^max]  denote  the  current  interval,  and  let  <p,  denote  the  interior  angle  Zp,_iPiP,+i 
for  0  <  t  <  n  —  1.  (For  convenience,  arithmetic  operations  performed  on  the  subscripts  of 
vertices  or  internal  angles  are  regarded  as  followed  by  a  ‘  mod  n’  operation;  for  example,  p_i 
is  identified  with  p„_i  and  Pn  with  po.)  In  the  algorithm,  $ieft  and  bright  keep  track  of  the 
angle  between  p/pj+i  and  popi  and  the  angle  between  PrPr+i  and  popi  respectively.  The 
algorithm  proceeds  as  follows: 

Algorithm  1 

Step  1  Start  at  the  pose  such  that  edge  poPi  aligns  with  the  right  side  of  C.  Locate  the 
vertex  pi  in  contact  with  the  left  side  of  C.  Set  $ieft  ♦-  I37=/+i(^  ~  '^i)i  bright  *-  0, 
^min  I  -  2>  I  and  r  4-  0, 

Step  2  The  current  orientation  interval  has  its  left  endpoint  at  d^\^.  Output  the  elliptic 
curve  (now  parameterized  by  6)  resulting  from  sliding  edge  pipr  in  cone  C. 

Next  determine  the  right  endpoint  ^»..v  of  the  current  interval.  Let  xj;  be  the  angle 
intersected  by  rays  pi-ipi  and  prPr-i’,  set  ^  4—  nil  if  they  do  not  intersect.  There  eire 
three  different  cases: 

Case  1  Ip  <  <f>  or  rp  =  nil.  [Advance  p,  clockwise  to  the  next  vertex.]  Set  bright  ♦— 

bright  "I"  ^  ^rj  ^max  *  bright  d"  ^  2’  and  r  4  r  1. 

Case  2  xp  >  <f>.  [Advance  p/.]  Set  $ieft  *—  $ieft  TT  (pi,  ^max  *  ^left  d"  2  2’  and 

/4-/-I. 

Case  S  [Advance  both  p/  and  pr-j  Set  $ieft  4—  $ieft  -f  tt  —  <p/,  bright  ^ 

bright  d-  TT  -  cpr,  d^ax  ♦“  ^left  d"  ^  d"  f ,  /  f  -  1,  and  T  4-  r  -  1. 

Output  the  current  interval  [^muu^max]-  Set  9nm  ^max  and  repeat  Step  2  until  r 
changes  from  1  to  0. 

The  number  of  intervals  produced  by  the  above  algorithm  cannot  exceed  2n,  because  each 
loop  of  Step  2  eidvances  either  pr  to  Pr-i,  or  p/  to  p/_i,  or  both  to  pr-i  and  p/_i  respectively, 
and  because  there  are  2n  vertices  in  total  (n  each  for  pi  and  pr)  to  advance  before  returning 
to  the  initial  incident  vertices  po  and  p,-. 

We  can  easily  apply  the  above  algorithm  for  the  general  case  in  which  the  vertex  of  cone 
C  is  at  an  arbitrary  point  (xq,  yo)  and  the  axis  of  the  cone  forms  an  angle  Oq  with  the  y  axis. 
Each  generated  interval  [a,  now  needs  to  be  right  shifted  to  [o-f  /^d-^oj-  Let  q  =  pr  and 
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^  =  Pi  if  po  is  above  edge  pipr,  and  let  q  =  pt  and  q'  =  pr  otherwise.  Then  the  corresponding 
locus  (z,  y)  of  po  is  determined  as,  assuming  pi  ^  Pr  Po, 


=  f-kPo|  •  cos  Ipipoq  +  \piPt\  •  kc  • 

\  sin  <p  J 


cos6> 


sin  <f) 

/■  I  .  .  I  .  sin(fT^o)\  .  « 

+  ^^l9Po|  •  sm  IpiPoq  -  \pipr\  •  A, - — j  sin^  +  zq; 


y  =  (^-\qpo\  •  sin  Ipipoq  ±  |p/Pr|  *  Ac 


cos(f  T^o) 


sin  (f> 


) 


cosd 


/|  I  ,  .  .  ,  cos(fT^o)\  .  « 

-  l^lqpol  •  cos  Ipipoq  ±  |p/pr|  •  A, - - jsin9  +  yo, 


where  Ac  =  cos(Zpipo9  +  Zpo??'  —  T  f )  ~  sin(^piPo9  +  ^Poqq'  -  T  f)-  Here  ‘±’ 

and  ‘t’  denote  ‘+’  and  or  ’  and  ‘+’,  according  as  po  is  above  or  below  edge  p/p^. 


2.2.3  Upper  Bounds 

The  preceding  subsection  tells  us  that  the  set  of  possible  poses  for  a  convex  polygon  in¬ 
scribed  in  one  cone  can  be  described  by  a  continuous  and  piecewise  elliptic  curve  defined 
over  orientation  space  [0,2jr).  We  call  this  curve  the  locus  curve  for  the  inscription.  This 
subsection  will  show  that  only  finite  possible  poses  exist  for  a  convex  n-gon  P  inscribed  in 
two  cones,  so  long  as  the  vertices  of  the  cones  do  not  coincide.  An  upper  bound  on  the 
number  of  possible  poses  can  be  obtained  straightforwardly  by  intersecting  two  locus  curves, 
each  resulting  from  the  inscription  of  P  in  one  cone. 

Claim  1  There  exist  no  more  than  8n  possible  poses  for  a  convex  n-gon  P  inscribed  in  two 
cones  Cl  and  C2  with  distinct  vertices. 

Proof.  Let  po, . . .  ,Pn-i  be  the  vertices  of  P  in  counterclockwise  order;  then  a  pose  of  P 
can  be  represented  by  the  location  of  po  as  well  as  the  angle  0  between  directed  edge  popl 
and  the  z  axis.  Let  Ci{d)  and  C2{$)  be  the  two  locus  curves,  for  the  inscriptions  of  P  in 
cones  Cl  and  C2  respectively.  We  only  need  to  show  that  Ci  and  C2  meet  at  most  8n  times, 
that  is,  they  pass  through  common  points  at  no  more  than  8n  values  of  9. 

It  is  known  that  each  C,-  consists  of  at  most  2n  elliptic  curves  defined  over  a  sequence  of 
intervals  that  partition  [0, 27r).  Intersecting  these  two  sequences  of  intervals  gives  a  partition 
that  consists  of  at  most  4n  intervals.  Within  each  interval  the  possible  orientations  (hence 
the  possible  poses)  of  P  can  be  found  by  computing  where  the  corresponding  pair  of  elliptic 
curves  meet.  According  to  the  last  subsection,  this  pair  of  curves  may  be  written  in  the 
parameterized  forms 


(a, a;  cos  0  -f  bix  sin  9  +  z,-,  a,j,  cos  9  +  biy  sin  9  -f  p,), 

for  i  =  1, 2.  Here  (z,-,  yi)  is  the  vertex  of  cone  C7,-,  and  a,>,  &,>,  a,„,  biy  are  constants  determined 
by  P  and  Ci.  Using  the  condition  (zi,pi)  ^  (2:2, 1/2),  we  suppose  zi  /  X2  without  loss  of 
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generality,  and  let  A  =  yf{aix  —  03*)^  +  (ii*  —  62*)^.  Then  it  is  not  hard  to  show  that  these 
two  curves  do  not  meet  if  |xi  —  X2I  >  A.  Otherwise  they  may  meet  only  at 

6  =  fi  —  a  and  9  =  ir  —  13  —  a, 

where  a =atan(  21*^22*.^  *^»~*»*)  and  □ 

The  upper  bound  8n  is  not  tight:  A  lower  one  can  be  obtained  even  without  using  two 
cones  to  constrain  the  polygon.  Notice  in  the  proof  above  that  the  bound  came  from  a 
partition  of  orientation  space  [0, 2r)  into  up  to  4n  intervals  which  combined  the  individual 
pose  constraints  imposed  by  the  two  cones.  Therefore  an  improvement  on  that  bound  must 
require  a  different  partitioning  of  [0, 2ir).  To  see  this,  we  regard  each  cone  as  the  intersec¬ 
tion  of  two  half-planes  and  decompose  its  constraint  on  the  polygon  into  two  constraints 
introduced  by  the  half-planes  independently. 

A  polygon  P  is  said  to  be  embedded  in  a  half-plane  A  if  P  is  contained  in  h  and  supported 
by  its  boimding  line.  Thus  P  is  inscribed  in  a  cone  if  and  only  if  it  is  embedded  in  the 
two  half-planes  defining  the  cone  by  intersection.  Two  cones  with  distinct  vertices  together 
provide  three  or  four  half-planes,  of  which  no  three  have  concurrent  bounding  lines,  i.e., 
bounding  lines  that  pass  through  a  common  point.  Such  three  half-planes  are  indeed  enough 
to  bound  the  number  of  possible  poses  of  P  within  6n. 

Theorem  1  There  exist  no  more  than  6n  possible  poses  for  a  convex  n-gon  P  embedded  in 
three  half-planes  with  non-concurrent  bounding  lines;  furthermore,  this  upper  bound  is  tight. 

Proof.  Let  /i,  I2  and  I3  be  the  bounding  lines  of  the  three  half  planes  respectively.  We  can 
assume  that  these  lines  are  not  all  parallel;  otherwise  it  is  easy  to  see  that  no  feasible  pose 
for  P  exists.  So  suppose  li  and  I2  intersect;  their  corresponding  half-planes  form  a  cone  in 
which  P  is  inscribed.  Let  the  orientation  of  P  be  represented  by  the  angle  9  between  the 
X  axis  and  some  directed  edge  of  P.  Then  orientation  space  [0,2ir)  is  partitioned  into  at 
most  2n  intervals,  according  to  which  pair  of  vertices  are  possibly  on  lx  and  I2  respectively. 
In  the  mean  time,  it  is  also  partitioned  into  exactly  n  intervals,  according  to  which  vertex 
is  possibly  on  I3.  Intersecting  the  intervals  in  these  two  partitions  yields  a  finer  partition  of 
[0,27r)  that  consists  of  at  most  3n  intervals,  each  containing  orientations  at  which  P  is  to  be 
supported  by  /i,  I2  and  I3  at  the  same  three  vertices. 

Let  us  look  at  one  such  interval,  and  let  p,,  pj  and  p*  be  the  vertices  of  P  on  /i,  I2  and  I3 
respectively  whenever  a  possible  orientation  exists  in  the  interval.  The  possible  orientations 
occur  exactly  where  I3  crosses  an  elliptic  curve  C(9)  traced  out  by  pk  when  sliding  p,-  and 
Pj  on  lx  and  I2  respectively.  Now  we  prove  that  this  interval  contains  at  most  two  possible 
orientations.  Note  C  does  not  degenerate  into  a  point  because  the  case  pi  =.  p^  —  pk  will 
never  happen,  given  lx,  I2  and  I3  are  not  concurrent.  Therefore  C  is  either  an  elliptic  segment 
monotonic  in  d  or  a  line  segment  that  attains  any  point  for  at  most  two  9  values.  (See 
Section  2.2.1.)  In  both  cases,  it  is  clear  that  C  crosses  I3  for  no  more  than  two  9  values.  Thus 
there  are  at  most  6n  possible  poses  in  orientation  space  [0,27r). 

Appendix  A  gives  an  example  in  which  a  polygon  can  actually  have  6n  poses  when 
embedded  in  three  given  half-planes,  thereby  proving  the  tightness  of  this  upper  bound.  □ 
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When  two  of  the  three  non-concurrent  bounding  lines  are  parallel,  we  can  similarly  derive 
a  smaller  but  tight  upper  bound  4n,  using  the  same  interval  partitioning  technique.  We  omit 
the  details  of  the  derivation  but  give  a  worst-case  example  in  Appendix  A. 

Since  2my  two  cones  with  distinct  vertices  are  formed  by  three  or  four  half-planes  with 
non-concurrent  bounding  lines,  amd  since  embedding  a  polygon  in  three  half-planes  with 
non-concurrent  bounding  lines  is  equivalent  to  inscribing  it  in  any  two  cones  determined  by 
intersecting  a  pair  of  the  half-planes,  we  immediately  have 

Corollary  1  There  exist  at  most  6n  possible  poses  for  a  convex  n-gon  inscribed  in  two  cones 
with  distinct  vertices,  and  this  upper  bound  is  tight. 

Would  more  half-planes  (or  cones)  further  reduce  the  number  of  possible  poses  for  an 
embedded  polygon  to  be  asymptotically  less  than  n?  The  answer  is  no.  For  example,  an 
embedded  regular  n-gon  will  always  have  kn  possible  poses,  where  1  <  <  6,  no  matter  how 

many  half-planes  are  present.  The  experimental  results  in  Section  2.3  will  show  that  two 
cones  (or  four  half-planes)  are  usually  sufficient  to  determine  a  unique  pose  for  the  polygon. 

2.2.4  An  Algorithm  for  Inscription 

With  the  results  in  the  previous  subsections,  we  here  present  an  algorithm  that  computes 
all  possible  poses  for  a  convex  n-gon  P  to  be  inscribed  in  m  cones  Ci, . . . ,  Cm,  where  m  >  2. 
(The  vertices  of  these  cones  are  assumed  to  be  distinct.)  Let  po,  •  •  •  ,Pn-i  he  the  vertices  of 
P  in  counterclockwise  order. 

Algorithm  2 

Step  1  [Compute  an  initial  set  of  poses  w.r.t.  two  cones.]  Solve  for  all  possible  poses  of  P 
when  inscribed  in  cones  Ci  and  C2  (use  Algorithm  1  and  see  the  proof  of  Claim  1), 
and  let  set  S  consist  of  the  resulting  poses  (edready  sorted  by  orientation).  Set  i  ♦—  3. 

Step  2  [Verify  with  the  remaining  cones.]  If  i  =  m  -f  1  or  5  =  0  then  terminate.  Otherwise 
go  to  Step  3  if  |5|  =  1  or  |5|  =  2.  Otherwise  apply  Algorithm  1  to  generate  the  locus 
curve  Ci(0)  for  the  inscription  of  P  in  C,-.  Sequentially  verify  whether  each  pose  in  S 
is  on  Ci{6),  deleting  from  5  those  poses  that  are  not.  Set  i  <—  t  -|- 1  and  repeat  Step  2. 

Step  3  [More  efficiently  verify  one  or  two  poses.]  For  each  pose  in  5,  let  polygon  P'  be  P 
in  that  pose  and  do  the  following:  For  i  <  j  <  m  construct  the  supporting  cone  Cj 
of  P'  at  the  vertex  of  cone  Cj;  if  there  exists  some  Cj  that  does  not  coincide  with  Cj, 
then  delete  the  corresponding  pose  from  S. 

When  the  above  algorithm  terminates,  set  S  will  contjiin  all  possible  poses  for  the  in¬ 
scription.  Corollary  1  shows  that  there  are  at  most  6n  poses  in  S  after  Step  1.  Since 
the  supporting  cone  of  P  from  a  point  can  be  constructed  in  time  O(logn)  using  binary 
search  [48],  the  running  time  of  the  algorithm  is  0{{k  —  l)n  -|-  (m  —  fc  -|- 1)  log  n),  i.e.,  0(mn) 
in  the  worst  case,  where  k  is  the  value  of  variable  i  when  leaving  Step  2.  However,  the 
experiments  in  Section  2.3  will  demonstrate  that  k  =  Z  almost  always  holds,  hence  Step  2 
will  almost  never  get  executed  more  than  once,  reducing  the  running  time  to  0{n  +  m  log  n). 
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Since  a  convex  m-gon  Q  is  naturally  the  intersection  of  m  cones,  each  with  a  vertex 
of  Q  as  its  vertex  and  with  the  internal  angle  at  that  vertex  as  its  apex  angle,  we  can 
use  Algorithm  2  to  compute  all  possible  inscriptions  of  P  in  Q  with  the  same  time  cost. 
This  problem  shall  be  called  the  polygon  inscription  problem,  which  can  be  regarded  as 
another  version  of  multi-cone  inscription  because  the  intersection  of  multiple  cones  is  always 
a  polygon  (possibly  unbounded  or  empty). 

2.3  Experiments 

The  first  set  of  experiments  were  conducted  to  find  out  how  many  possible  poses  usually 
exist  for  a  polygon  embedded  in  three  half-planes  with  non-concurrent  bounding  lines.  The 
results  are  summarized  in  Table  1. 


#  tests 

data 

source 

#  poly  vertices 

#  possible  poses 

range 

mean 

range 

mean 

10  sq. 

3-9 

5.9376 

1-16 

5.5436 

10  cir. 

6.1208 

2-16 

5.7006 

100  sq. 

6-19 

11.917 

2-16 

7.242 

1000 

100  cir. 

10-22 

15.108 

2-24 

10.012 

1000  sq. 

11-24 

18.43 

2-18 

7.84 

■ciiiiiggni 

23-43 

33.595 

2-60 

18.709 

1000 

cir.  mar. 

3-15 

8.98 

2-14 

4.537 

Table  1:  Experiments  on  embedding  a  polygon  in  three  half-planes. 

Seven  groups  of  convex  polygons  were  tested.  The  first  six  groups  consisted  of  convex 
hulls  generated  over  10,  100  and  1000  random  points  successively,  and  for  each  number,  in 
two  kinds  of  uniform  distributions — inside  a  square  and  inside  a  circle  respectively.  It  can 
be  seen  in  the  table  that  the  polygons  in  these  groups  had  a  wide  range  (3-43)  of  sizes, 
i.e.,  numbers  of  vertices,  but  their  shapes  were  not  arbitrary  enough,  approaching  either  a 
square  or  a  circle  when  large  number  of  random  points  were  used.  So  we  introduced  the  last 
group  of  data  that  consisted  of  polygons  generated  by  a  method  called  circular  march  which 
outputs  the  vertices  of  a  convex  polygon  as  random  points  inside  a  circle  in  counterclockwise 
order.  The  size  of  a  polygon  in  this  group  was  randomly  chosen  between  3  and  15. 

Given  a  convex  polygon,  three  supporting  lines,  each  bounding  a  half-plane  on  the  side  of 
the  polygon,  were  generated  according  to  the  uniform  distribution;  namely,  with  probability 
each  line  passed  through  vertex  p,-  with  internal  angle  <p,.  An  additional  check  was 
performed  to  ensure  that  these  lines  were  not  attached  to  the  same  vertex  of  the  polygon. 
The  number  of  possible  poses  for  the  polygon  to  be  embedded  in  these  generated  half-planes 
was  then  computed,  and  the  summarized  results  for  all  group  are  listed  in  the  last  two 
columns  of  the  table. 

Table  1  tells  us  that  three  half-planes  are  insufficient  to  limit  all  possible  poses  of  an 
embedded  polygon  to  a  unique  one,  namely,  the  real  pose;  in  fact  the  table  suggests  that 
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Figure  5:  flight  possible  poses  for  a  convex  6-gon  bounded  by  three  supporting  lines  (as  taken 
from  a  sample  run).  The  first  one  represents  the  real  pose  whose  supporting  lines  as  shown  were 
generated  randomly;  the  remaining  seven  represent  all  other  poses  consistent  with  the  supporting 
line  constrmnts.  Notice  in  this  example  that  three  of  the  eight  poses  occurred  when  a  vertex  of  the 
pdiygon  coincided  with  an  intersection  of  two  supporting  lines. 

linear  (in  the  size  of  the  polygon)  number  of  possible  poses  will  usually  exist.  We  can  see 
in  the  table  that  despite  the  appearances  of  cases  with  one  or  two  possible  poses,  the  ratio 
between  the  mean  of  numbers  of  possible  poses  £md  mean  polygon  size  lies  in  the  approximate 
range  0.43-0.93,  decreasing  very  slowly  as  the  mean  polygon  size  in  a  group  increases.  These 
results  tend  to  support  a  conjecture  that  in  the  average  case  there  exist  0{n)  possible  poses 
for  a  convex  n-gon  embedded  in  three  half-planes  with  non-concurrent  bounding  lines. 

The  above  conjecture  may  be  very  difficult  to  prove.  However,  a  plausible  explanation 
for  the  experimental  results  can  be  sought.  Recall,  a  polygon  with  three  half-planes  defines 
a  partition  of  orientation  space  [0,27r)  into  at  most  3n  intervals,  each  conteuning  orientations 
that  would  allow  the  same  three  incident  vertices  whenever  a  possible  pose  exists  at  that 
orientation.  The  feasible  orientations  in  each  interval  occur  when  one  supporting  line  crosses 
the  locus  curve  of  its  associated  incident  vertex.  The  locus  curve  results  from  moving  the 
other  two  incident  vertices  along  their  supporting  lines.  As  these  curves  (for  all  intervals) 
may  often  cluster  together,  the  likelihood  that  they  get  crossed  0{n)  times  in  total  by  the 
first  supporting  line  is  quite  large.  This  happened  particularly  often  during  the  experiments 
when  a  vertex  coincided  with  an  intersections  of  two  supporting  lines.  (See  Figure  5.) 

The  purpose  of  the  second  set  of  experiments  was  to  study  how  many  poses  usually  exist 
for  a  convex  polygon  inscribed  into  two  or  more  cones  with  distinct  vertices.  We  first  tested 
with  two  cones  using  the  same  source  of  random  data  generated  in  the  way  we  did  for  the 
first  set  of  experiments,  and  the  results  are  shown  in  Table  2.  Since  a  polygon  w£is  always 
generated  inside  a  square  (circle),  the  cone  vertices  were  chosen  aus  random  points  uniformly 


#  tests 

data 

source 

#  poly  vertices 

^  possible  poses 

ISl 

mean 

range 

mean 

10000 

■ligf 

ESOi 

1-2 

1.036 

10000 

3-10 

6.1113 

1-2 

1.0117 

1000 

100  sq. 

6-18 

12.002 

1-2 

1.01 

1000 

100  cir. 

9-21 

15.138 

1-1 

1 

1000 

10-26 

18.073 

1-2 

1.002 

1000 

26-45 

33.665 

1-2 

1.003 

1000 

3-15 

9.009 

1-2 

1.106 

Table  2:  Ebcperiments  on  inscribing  polygons  with  two  cones. 

distributed  between  this  boundary  and  a  larger  square  (circle).  The  ratio  between  these  two 
squares  (circles)  was  set  uniformly  to  be  |  for  all  seven  groups  of  data. 

In  contrast  to  Table  1,  Table  2  tells  us  that  two  cones  allow  a  unique  pose  of  an  inscribed 
polygon  in  most  cases.  In  each  group  of  tests,  only  cases  with  one  pose  or  two  poses  occurred, 
and  the  mean  of  possible  poses  stayed  very  close  to  1,  independent  of  the  mean  polygon  size. 
(It  is  not  hard  for  us  to  see  that  the  percentage  of  two-pose  cases  was  very  low  in  the  range 
0%-3.6%  for  the  first  six  groups  of  data.  The  pecentage  10.6%  for  the  seventh  group  was  a 
bit  high  but  expected,  because  polygons  generated  by  circular  march  were  more  likely  to  be 
in  a  certain  shape  that  would  often  incur  two  pc»sible  poses,  as  we  will  discuss  later.) 

Tests  were  also  conducted  with  3-10  cones  on  reproduced  data  of  four  of  the  seven  groups, 
while  the  other  experiment  parameters  were  kept  the  same.  As  shown  in  Table  3,  the  means 
of  possible  poses  did  not  decrease  dramatically  as  compared  to  those  in  Table  2.  Finally, 


#  tests 

data 

source 

#  poly  vertices 

#  possible  poses 

mean 

range 

mean 

10000 

10  sq. 

3-10 

5.95 

1-2 

1.0056 

1000 

100  sq. 

■ani 

11.901 

1-1 

1 

10-27 

18.26 

1-1 

1 

1000 

cir.  mar. 

3-15 

8.96 

1-2 

1.022 

10000 

Msm 

5.9587 

1-2 

10  sq. 

3-10 

5.9741 

1-2 

1.1738 

Table  3:  More  experiments  on  inscribing  polygons  with  cones.  The  first  four  groups  of  data 
were  tested  with  a  random  number  (between  3  and  10)  of  cones;  the  last  two  groups  were  tested 
with  two  cones  with  vertices  chosen  from  wider  ranges — the  ratio  between  the  two  squares  (circles) 
determining  the  range  of  cone  vertices  was  set  to  |  and  ^  respectively  for  these  two  groups. 

we  repeated  the  first  group  of  tests  with  two  cones  but  chose  their  vertices  from  two  wider 
ranges  (with  the  previous  ratio  |  replaced  by  ratios  j  and  ^  respectively),  and  the  results 
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Figure  6:  Two  possible  poses  for  a  convex  6-gon  inscribed  in  two  cones  (as  taken  from  a  sample 
run).  The  two  cones  are  supporting  the  polygon  at  a  pair  of  vertices. 

are  shown  in  the  last  two  rows  of  the  same  table. 

The  experimental  result  that  two  non-incident  cones  usually  allow  a  unique  pose  of  an 
inscribed  convex  n-gon  P  has  in  fact  a  very  intuitive  explanation.  As  mentioned  before,  two 
such  cones  generally  provide  four  half-planes,  any  three  of  which  will  limit  the  number  of 
possible  poses  of  P  to  at  most  6n.  Let  polygons  where  m  <  6n,  represent  P 

in  all  possible  poses  respectively  when  embedded  in  the  first  three  hjilf-planes;  then  those 
Pi  corresponding  to  the  final  possible  poses  must  be  supported  by  the  bounding  line  /  of 
the  fourth  half-plane.  So  the  probability  that  a  two-pose  case  occurs  is  no  more  than  the 
probability  that  /  passes  through  a  vertex  of  Pi  and  a  vertex  of  Pj,  for  i  ^  j.  Note  the 
vertices  of  Pi,..., Pm  together  occupy  ©(mn)  points  in  the  plane  in  general,  only  one  of 
which  must  lie  on  /.  If  no  two  of  these  vertices  coincide,  the  probability  that  I  passes 
through  two  vertices  of  different  polygons  is  zero  (assuming  that  /  is  independent  of  the 
other  three  bounding  lines),  which  means  that  a  two-pose  case  almost  never  occurs  in  this 
situation.  Otherwise  suppose  two  vertices  of  Pi  and  Pj  respectively  are  at  the  same  point 
p  for  some  i  j,  then  the  probability  that  /  passes  through  p  is  0(;^).  This  is  Q(^)  in 
the  average  case,  as  suggested  that  m  =  0(n)  by  the  first  set  of  experiments.  Since  in  the 
usual  case  there  only  exist  a  constant  number  of  such  coincident  vertex  pairs,  the  probability 
is  an  approximate  upper  bound  on  how  often  two-pose  cases  occur.  This  bound  turns 
out  to  be  consistent  with  the  percentages  of  two-pose  cases  in  Table  2. 

It  was  observed  during  the  experiments  that  a  large  number  of  two-pose  cases  occurred 
when  both  cones  happened  to  be  supporting  the  polygon  at  the  same  pair  of  vertices.  (See 
Figure  6.)  The  two  possible  orientations  differed  by  ir,  and  each  supporting  vertex  in  one  pose 
coincided  with  the  other  in  the  other  pose.  This  situation  often  happened  when  the  distance 
between  one  pair  of  vertices  of  the  polygon  was  much  larger  than  the  distance  between  any 
other  pair  of  vertices,  or  when  the  sites  were  far  away  from  the  polygon  (as  evidenced  by  the 
high  percentages  of  two-pose  cases  in  the  last  two  groups  of  tests  in  Table  3). 


3  Sensing  by  Point  Sampling 

3.1  The  Point  Sampling  Problem 

Often  the  possible  poses  in  which  a  part  settles  on  the  assembly  table  axe  not  continuous  but 
of  a  small  number,  either  reduced  by  a  parts  feeder  or  limited  by  the  mechanical  constraints 
imposed  by  a  sequence  of  planned  manipulations.  (See  Erdmann  and  Mason’s  tray-tilting 
method  [22]  and  Brost’s  squeeze-grasping  method  [10]  for  examples  of  the  latter  case.)  These 
facts  together  allow  the  implementation  of  more  effective  sensing  mechanisms.  The  efficiency 
of  such  sensing  mechanisms  depends  on  both  the  time  cost  of  the  physical  operations  and 
the  time  complexity  of  the  algorithms  involved.  Consequently,  minimizing  one  or  both  of 
these  two  factors  has  become  an  important  aspect  of  sensor  design. 

In  order  to  illustrate  the  goals  of  this  section,  consider  a  polygonaJ  part  resting  on  a 
horizontal  assembly  table.  The  table  is  bounded  by  vertical  fences  at  its  bottom-left  corner, 
as  shown  in  Figure  7.  Pushing  the  part  towards  that  corner  will  eventually  cause  the  part 
to  settle  in  one  of  the  12  stable  poses  listed  in  the  figure.  (Note  to  reach  a  stable  pose  both 
fences  must  be  in  contact  with  some  vertices  of  the  part  while  at  least  one  fence  must  be 
in  contact  with  no  less  than  two  vertices.)  In  order  to  distinguish  between  these  12  poses, 
the  robot  has  marked  4  points  on  the  table  beforehand,  so  it  can  infer  the  pose  from  which 
marks  are  covered  by  the  part  and  which  are  not*. 

The  above  “shape  recovery”  method  is  named  sensing  by  point  sampling,  as  a  loose 
analogy  to  the  reconstruction  of  band  limited  functions  by  sampling  on  a  regular  grid  in 
signal  processing.  To  save  the  expense  of  sampling,  the  robot  wants  to  mark  as  few  points 
as  possible.  The  problem:  How  to  compute  a  minimum  set  of  points  to  be  marked  so  that 
parts  of  different  types  and  poses  can  be  distinguished  from  each  other  by  this  method? 

3.1.1  Related  Work 

Orienting  mechanical  parts  was  early  studied  by  Grossman  and  Blasgen  [31].  They  used  a 
vibrating  box  to  constrain  a  part  to  a  small  finite  number  of  possible  stable  poses,  cind  then 
determined  the  particular  pose  by  a  sequence  of  probes  using  a  tactile  sensor.  Natarajan  [46] 
examined  a  similar  strategy  of  detecting  the  orientations  of  polygonfil  and  polyhedral  objects 
with  an  analysis  of  the  numbers  of  sensors  sufficient  and  necessary  for  the  task. 

More  recent  related  work  includes  [8]  and  [3].  [8]  shows  that  the  problem  of  deciding 
whether  k  line  probes  are  sufficient  to  distinguish  a  convex  polygon  from  a  collection  of  n 
convex  polygons  is  NP-complete.  This  result  is  very  similar  to  our  Theorem  3.  A  variation  of 
the  line-probing  result  in  [8]  would  give  us  the  point  sampling  result  of  Theorem  3.  [3]  proves 
a  similar  result  as  well,  namely  that  the  problem  of  constructing  a  decision  tree  of  minimum 
height  to  distinguish  among  n  polygons  using  point  probes  is  NP-complete.  This  result  holds 
even  if  all  the  polygons  are  convex.  [3]  also  exhibits  a  greedy  approximation  algorithm  for 
constructing  such  a  decision  tree.  This  result  is  simileur  to  our  approximation  algorithm  of 
Section  3.4,  with  a  similar  ratio  bound.  The  difference  is  that  our  greedy  algorithm  seeks  to 

^This  can  be  implemented  in  multiple  ways,  such  as  placing  light  detectors  in  the  table,  probing  at  the 
points,  or  if  the  robot  has  a  vision  system,  taking  a  scene  image  and  checking  the  corresponding  pixel  values. 
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(a) 


(b) 

Figure  7:  Sensing  by  point  sampling,  (a)  The  12  possible  stable  poses  of  an  assembly  part 
after  pushing,  along  with  4  sampling  points  (optimal  by  Lemma  1)  found  by  our  implementation 
of  the  approximation  algorithm  to  recognize  these  poses,  where  dashed  line  segments  are  fences 
perpendicular  to  the  plane,  (b)  The  planar  subdivision  formed  by  these  poses  which  consists  of  610 
regions. 
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minimize  the  total  number  of  probe  points  rather  than  the  tree  height.^ 

Closely  related  work  includes  the  research  by  Romanik  and  others  on  geometric  testability 
(see,  for  example,  [50],  [51],  and  [1]).  Their  research  develops  strategies  for  verifying  a  given 
polygon  using  a  series  of  point  probes.  Moreover,  the  research  examines  the  testability  of 
more  general  geometric  objects,  such  as  polyhedra,  and  develops  conditions  that  determine 
whether  a  class  of  objects  is  (approximately)  testable. 

A  number  of  researchers  have  looked  at  the  problem  of  determining  or  distinguishing 
objects  using  finger  probes.  Finger  probing  is  closely  related  to  sensing  by  point  sampling, 
as  indicated  by  our  discussion  of  [8].  For  a  more  extensive  survey  of  probing  problems  and 
solutions  see  the  paper  by  Skiena  [54]. 

There  would  seem  to  be  connections  between  our  work  and  the  concept  of  VC-dimension 
often  used  in  learning  theory.  For  instance,  in  this  section  we  develop  the  notion  of  a 
“discriminating  set”  to  distinguish  different  polygons.  The  concept  of  a  discriminating  set 
bears  some  resemblance  to  the  idea  of  shattered  sets  associated  with  VC-dimension.  However, 
discriminating  sets  and  shattered  sets  are  different.  A  minimum  discriminating  set  is  the 
smallest  set  of  points  that  uniquely  identifies  every  object  in  a  set  of  objects,  whereas  VC- 
dimension  is  the  size  of  the  largest  set  of  points  shattered  by  the  set  of  objects.  Thus,  the 
yC-dimension  of  a  finite  class  gives  a  lower  bound  on  the  size  of  a  minimum  discriminating 
set.  For  dense  polygon  distributions,  the  two  cardinalities  will  be  the  same,  namely  logn, 
where  n  is  the  number  of  polygons.  For  sparsely  distributed  polygons,  the  two  cajdiniilities 
are  different.  For  instsmce,  the  VC-dimension  can  be  1,  while  the  minimum  discriminating 
set  has  size  n  —  1.  See  Figure  8. 

Finally,  the  work  described  in  this  section  is  part  of  our  larger  research  goal  to  un¬ 
derstand  the  information  requirements  of  robot  tasks.  Related  work  includes  the  sensor 
design  methodology  of  Erdmann  [21]  and  the  information  invariants  of  Donald  et  al.  [19]. 
[21]  proposes  a  method  for  designing  sensors  based  on  the  particuleur  manipulation  task  at 
hand.  The  resulting  sensors  satisfy  a  minimality  property  with  respect  to  the  given  task 
goal  and  the  available  robot  actions.  [19]  investigates  the  relationship  between  sensing,  ac¬ 
tion,  distributed  resources,  communication  paths,  and  computation,  in  the  solution  of  robot 
tasks.  That  work  provides  a  method  for  comparing  disparate  sensing  strategies,  and  thus 
for  developing  minimal  or  redundant  strategies,  as  desired. 

3.1.2  The  Formal  Problem 

Consider  n  simple  polygons  in  the  plane,  not  necessarily  disjoint  from  each  other. 

We  wish  to  locate  the  minimum  number  of  points  in  the  plaine  such  that  no  two  polygons  Pi 
and  Pj,  j,  contain  exactly  the  same  points.  In  order  to  avoid  ambiguities  in  sensing,  we 
require  that  none  of  the  located  points  lie  on  any  edge  of  Pi, ... ,  P„.  The  planar  subdivision 
formed  by  Pi, . . . ,P„  divides  the  plane  into  one  unbounded  region,  some  bounded  regions 

^It  is  easy  to  give  an  example  for  which  a  minimum  height  decision  tree  uses  more  than  minimum  number 
of  total  probes,  while  a  decision  tree  with  minimum  number  of  total  probes  does  not  attain  the  minimum 
height.  Consider  the  problem  of  discriminating  sets  {a,6,a',e'},  {a, a'},  {6,6',d'},  {c,6'},  {d,d}  and  0  which 
can  be  viewed  as  probing  a  collection  of  polygons  by  the  later  transformation  technique  in  the  proof  of 
Theorem  3.  The  decision  tree  using  minimum  probes  a,b,c,d  always  has  height  4,  and  the  decision  tree 
using  probes  a\b',cf,d',e'  can  achieve  minimum  height  3- 
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outside  Pi, . . . ,  P„,  called  the  “holes”,  and  some  bounded  regions  inside.  (For  example,  the 
12  polygons  in  Figure  7(a)  form  the  subdivision  in  Figure  7(b)  which  consists  of  610  regions, 
none  of  which  is  a  hole.)  Immediately  we  make  two  observations:  (1)  Points  on  the  edges 
of  the  subdivision  or  in  the  interior  of  the  unbounded  region  or  in  a  “hole”  do  not  need  to 
be  considered  as  locations;  (2)  for  each  bounded  (open)  region  inside  some  polygon  only  one 
point  needs  to  be  considered. 

Let  n  denote  the  set  of  bounded  regions  in  the  subdivision  which  are  contained  in  at 
least  one  of  Pi, . . .  ,Pn.  Each  polygon  Pj,  1  <  i  <  n,  is  partitioned  into  one  or  more  such 
regions;  we  write  u;  C  p  when  a  region  u  is  contained  in  polygon  P,.  A  region  basis  for 
polygons  Pi , . . . ,  P„  is  a  subset  A  C  such  that 

{  u;  I  u;  €  A  and  ci;Cp}^{c<;|u;€A  and  a;  C  Pj  }, 

for  1  <  i  ^  j  <  n;  that  is,  each  P  contains  a  distinct  collection  of  regions  from  A.  A  region 
basis  A*  of  minimum  cardinality  is  called  a  minimum  region  basis.  Thus  the  problem  of 
sensing  by  point  sampling  becomes  the  problem  of  finding  a  minimum  region  basis  A*.  We 
will  call  this  problem  Region  Basis  and  focus  on  it  throughout  the  section.  The  following 
lemma  gives  the  upper  and  lowe..'  bounds  for  the  size  of  such  A*. 

Lemma  1  A  minimum  region  basis  A*  forn  polygons  Pi,...,P„  satisfies  [logn]  <  |A*|  < 
n  —  1. 

Proof.  To  verify  the  lower  bound  flog  n] ,  note  that  each  of  the  n  polygons  must  contain 
a  distinct  subset  of  A*;  so  n  <  the  cardinality  of  the  power  set  2^*. 

To  verify  the  upper  bound  n  —  1,  we  incrementally  construct  a  region  basis  A  of  size  at 
most  n  —  1.  This  construction  is  similar  to  Natarajan’s  Algorithm  2  [46j.  Initially,  A  =  0.  If 
n  >  1,  without  loss  of  generality,  assume  Pi  hcis  the  smallest  area.  Then  there  exists  some 
region  wi  €  outside  Pi.  Split  {Pi, . . . ,  P„}  into  two  nonempty  subsets,  one  including  those 
p  containing  wi  and  the  other  including  those  not;  and  add  u;i  into  A.  Recursively  split 
the  resulting  subsets  in  the  same  way,  and  at  each  split,  add  into  A  its  defining  region  (as 
we  did  with  wi)  if  this  region  is  not  already  in  A,  until  every  subset  eventually  becomes  a 
singleton.  The  A  thus  formed  is  a  region  basis.  Since  there  are  n  —  1  splits  in  total  and  each 
split  adds  at  most  one  region  into  A,  we  have  |A|  <  n  —  1.  □ 

Figure  8  gives  two  examples  for  which  |A*|  =  flogn]  and  |A*j  =  n  —  1  respectively. 
Therefore  these  two  bounds  are  tight. 

We  can  view  all  the  bounded  non-hole  regions  as  elements  of  fl,  and  all  the  polygons 
Pi, . . .  ,P„  as  subsets  of  Q.  Then  a  region  basis  A  is  a  subset  of  fi  that  can  discriminate 
subsets  Pi,...,P„  by  intersection.  Hence  the  Region  Basis  problem  can  be  rephrased  as: 
Find  a  subset  of  H  of  minimum  size  whose  intersections  with  any  two  subsets  P,  and  Pj, 
1  ^  ^  i  ^  are  not  equal.  The  general  version  of  this  set-theoretic  problem,  in  which 

stands  for  an  arbitrary  finite  set  and  Pi, . . . ,P„  stand  for  arbitrary  subsets  of  fl,  we  call 
Discriminating  Set.  We  have  thus  reduced  Region  Basis  to  Discriminating  Set,  and  the 
former  problem  will  be  solved  once  we  solve  the  latter  one. 

Let  us  analyze  the  amount  of  computation  required  for  the  geometric  preprocessing  to 
reduce  Region  Basis  to  Discriminating  Set.  Let  m  be  the  total  size  of  Pi, ... ,  P„,  i.e.,  the  sum 
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(»)  (b) 

Figure  8:  Two  examples  whose  minimum  region  basis  sizes  achieve  the  lower  bound  flog  n]  and  the 
upper  bound  n—  1  respectively.  Bounded  regions  in  the  examples  are  labelled  with  numbers,  (a)  For 
1  <  t  <  n  polygon  Pi  is  defined  to  be  the  boundary  of  the  union  of  regions  flog  n1  + 1, . . . ,  flog  n"}  +t, 
and  all  regions  k  with  1  <  A;  <  flog  nl  such  that  the  A;th  bit  of  the  binary  representation  (radix  2) 
for  t  -  1  is  1.  Thus  A*  =  {1,2,...,  flog n]}.  (b)  The  polygons  Pi,...,P„  contain  each  other  in 
increasing  order:  A*  =  {2, 3, . . . ,  n}. 

of  the  number  of  vertices  each  polygon  has;  trivially  m  >  3.  Then  the  planar  subdivision 
these  polygons  define  has  at  most  s  vertices,  where  3  <  s  <  By  Euler’s  relation 

on  planar  graphs,  the  number  of  regions  and  the  number  of  edges  are  upper  bounded  by 
2s  —  4  and  3s  —  6  respectively.  So  we  can  construct  the  planar  subdivision  either  in  time 
0(m  log  m  +  s)  using  an  optimal  algorithm  for  intersecting  line  segments  by  Chazelle  and 
Edelsbrunner  [13],  or  in  time  C?(slogm)  using  a  simpler  plane  sweep  version  by  Nievergelt 
and  Preparata  [47].  To  obtain  the  set  of  regions  each  polygon  contains,  we  only  need  to 
traverse  the  portion  of  the  subdivision  bounded  by  that  polygon,  which  takes  time  0(s).  It 
follows  that  the  reduction  to  Discriminating  Set  can  be  done  in  time  0(m  log  m  +  ns),  or 
O(nm^)  in  the  worst  case. 

Here  is  a  short  summary  of  the  structure  of  Section  3:  Section  3.2  proves  the  NP- 
completeness  of  Discriminating  Set;  based  on  this  result.  Section  3.3  establishes  an  equiva¬ 
lence  between  Discriminating  Set  and  Region  Basis,  hence  proving  the  latter  problem  NP- 
complete;  Section  3.4  presents  an  0{n^m^)  approximation  algorithm  for  Region  Basis  with 
ratio  2 Inn  and  shows  that  further  improvements  on  this  ratio  are  hard;  and  Section  3.5 
closes  up  with  a  complexity  analysis  of  various  subproblems  of  Discriminating  Set,  along 
with  the  definition  of  a  family  of  related  NP-complete  problems  called  A:-Independent  Sets. 
We  have  implemented  our  approximation  algorithm  and  have  tested  it  on  both  real  data 
taken  from  mechanical  parts  and  random  data  extracted  from  the  arrangements  of  random 
lines.  The  algorithm  works  very  well  in  practice. 

3.2  Discriminating  Set 

Given  a  collection  C  of  subsets  of  a  finite  set  X,  suppose  we  want  to  identify  these  subsets  just 
from  their  intersections  with  some  subset  D  Q  X.  Thus  D  must  have  distinct  intersection 
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with  every  member  of  C,  that  is, 

D  r\  S  ^  D  r\T,  for  all  S,T  €.  C  and  S  ^  T. 

We  call  such  a  subset  D  a  discriminating  set  for  C  with  respect  to  X.  From  a  different  point 
of  view,  each  element  x  E  D  can  be  regarded  as  a  bineiry  “bit”  that,  to  represent  any  subset 
sex,  gives  value  ‘1’  if  x  €  5  and  value  ‘0’  otherwise.  In  such  a  way  D  represents  an 
encoding  scheme  for  subsets  in  C. 

Below  we  show  that  the  problem  of  finding  a  minimum  discriminating  set  is  NP-complete. 
As  usual,  we  consider  the  decision  version  of  this  minimization  problem; 

Discriminating  Set  (D-Set) 

Let  C  be  a  collection  of  subsets  of  a  finite  set  X  and  /  <  1X|  a  non-negative  integer.  Is  there 
a  discriminating  set  D  C  X  for  C  such  that  |D|  <  /? 

Our  proof  of  the  NP-completeness  for  D-Set  uses  a  reduction  from  Vertex  Cover  (VC)  which 
determines  if  a  graph  G  =  {V,E)  has  a  cover  of  size  not  exceeding  some  integer  /  >  0,  i.e.,  a 
subset  V  QV  that,  for  each  edge  (it,  v)  €  B,  contains  either  u  or  v.  The  reduction  is  based 
on  a  key  observation,  that  for  any  three  finite  sets  Si,  S2  and  S3, 

Si  n  S2  Si  n  S3  if  and  only  if  5i  n  (52  A  S3)  ^  0, 

where  ‘A’  denotes  the  operation  of  symmetric  difference,  i.e.,  52  A  5$  =  (52  \  S3)  U  (53  \  52). 

Theorem  2  Discriminating  Set  is  NP-complete. 

Proof.  That  D-Set  €  NP  is  trivial. 

Next  we  establish  VC  ocp  D-Set,  that  is,  there  exists  a  polynomial-time  reduction  from  VC 
to  D-Set.  Let  G  =  (V,  E)  and  integer  0  <  /  <  |  V|  be  an  instance  of  VC.  We  need  to  construct 
a  D-Set  instance  (X,  C)  such  that  the  collection  C  of  subsets  of  X  has  a  discriminating  set 
of  size  /'  or  less  if  and  only  if  G  has  a  vertex  cover  of  size  I  or  less. 

The  construction  uses  the  component  design  technique  described  by  Garey  and  John¬ 
son  [25].  It’s  rather  natural  for  us  to  begin  by  including  every  vertex  of  G  in  set  X,  and 
assigning  each  edge  e  =  (u,v)  a  subset  5(e)  in  C  which  contains  at  least  u  and  u;  in  other 
words,  we  have  V  CX  and 

{it,u}  C  5(e)  €  C,  for  all  e  =  (u,i;)  6  E. 

In  order  to  ensure  that  any  discriminating  set  D  for  C  contains  at  least  one  of  u  and  v  from 
subset  5(e),  we  add  an  auxiliary  subset  into  C  that  consists  of  some  new  elements  not  in 
V,  and  in  the  meantime  define 

5(e)  =  {u,v}  U  Ae. 

Hence  5(e)  A  Ae  =  {u,u};  and  D  D  {u,u}  ^  0  follows  directly  from  D  D  5(e)  ^  D  D  Ae. 
Since  any  discriminating  set  D'  for  {  Ae  |  e  €  }  can  also  distinguish  between  5(ei)  and 

5(c2),  and  between  5(ei)  and  Ae^,  for  any  ei,e2  €  E  and  ei  /  e2,  D'  unioned  with  a  vertex 
cover  for  G  becomes  a  discriminating  set  for  C.  Conversely,  every  discriminating  set  D  for 
C  can  be  split  into  a  discriminating  set  for  {  Ae  |  e  €  }  and  a  vertex  cover  for  G. 
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The  m  =  |£|  auxiliary  subsets  should  be  constructed  in  a  way  such  that  we  can  easily 
determine  the  size  of  their  minimum  discriminating  sets  in  order  to  set  up  the  entire  D-Set 
instance.  There  is  a  simple  way:  We  introduce  m  elements  ui,  02, . . . ,  ^  V  into  X,  and 

define  subsets  Ae,  for  e€  E,  to  be 

{^2}}  •••»  {®m}> 

where  the  order  of  mapping  does  not  matter.  It’s  clear  that  there  are  m  minimum  discrimi¬ 
nating  sets  for  the  above  subsets:  {oi, . . .  ,0^}  \  {^i},  1  <  i  <  m. 

Setting  r  =  /  -h  m  —  1,  we  have  completed  our  construction  of  the  D-Set  instance  as 

X  —  VU{flij...,  (lX>  ■  ■  ■  )  ^ 

C  =  {5(c)|c€^}U{Ae|e€£;}. 

The  construction  can  be  carried  out  in  time  0{\V\  -h  |£|).  We  omit  the  remaining  task  of 
verifying  that  G  has  a  vertex  cover  of  size  at  most  I  if  and  only  if  C  has  a  discriminating  set 
of  size  at  most  /  -f-  m  —  1.  □ 

One  thing  about  this  proof  is  worthy  of  note.  All  subsets  in  C  constructed  above  have 
at  most  three  elements.  This  reveals  that  D-Set  is  still  NP-complete  even  if  IS]  <  3  for  all 
5  €  C,  a  stronger  assertion  than  Theorem  2.  The  subproblem  where  all  5  €  C  have  |5|  <  1 
is  obviously  in  P,  for  an  algorithm  can  simply  count  |C|  in  linear  time  and  then  answer  “yes” 
if  /  >  ICI  —  1  and  “no”  if  0  <  /  <  |C|  —  1.  For  the  remauning  case  in  which  all  5  6  C  have 
|5|  <  2,  we  will  prove  in  Section  3.5  that  the  NP-completeness  still  holds.  However,  the 
proof  will  be  a  bit  more  involved  than  the  one  we  just  gave  under  no  restriction  on  |5|. 

At  the  end  of  Section  3.2,  we  give  a  problem  that  is  equivalent  to  D*Set: 

Rx)w- Differing  Submatrix 

Given  an  m  x  n  matrix  A  of  O’s  and  I’s  and  integer  0  <  /  <  n,  is  there  an  m  x  /  matrix  B 
formed  by  /  columns  of  matrix  A  such  that  no  two  rows  of  B  are  identical? 


3.3  Region  Basis 

Now  that  we  have  shown  the  NP-completeness  of  D-Set,  the  minimum  region  basis  cannot 
be  computed  in  polynomial  time  through  the  use  of  an  efficient  algorithm  for  D-Set,  because 
no  such  algorithm  would  exist  unless  P  =  NP.  This  conclusion,  nevertheless,  leads  us  to 
conjecture  that  the  minimization  problem  Region  Basis  is  also  NP-complete.  Again  we 
consider  the  decision  version: 

Region  Basis  (RB) 

Given  n  polygons  Pi,..,,  Pn  and  integer  0  <  /  <  n  —  1,  does  there  exist  a  region  basis  A  for 
the  plansu:  subdivision  Cl  formed  by  Pi, . . . , such  that  | A|  <  /? 

The  condition  0  <  /  <  n  —  1  above  is  necessary  because  we  already  know  from  Lemma  1 
that  a  minimum  region  basis  has  size  at  most  n  —  1. 

Consider  a  mapping  P  from  the  set  of  RB  instances  to  the  set  of  D-Set  instances  that 
maps  regions  to  elements  and  polygons  to  subsets  in  a  one-to-one  manner.  Every  RB  instance 
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is  thus  mapped  into  an  equivalent  D-Set  instance,  as  pointed  out  in  Section  3.1.  We  claim 
that  T  is  not  onto.  Suppose  T  were  onto.  Then  the  elements  of  each  subset  in  a  D-Set 
instance  must  correspond  to  regions  in  some  RB  instance.  The  union  of  these  regions  must 
be  a  polygon,  and  this  polygon  must  map  to  the  subset  given  in  the  D-Set  instance.  However, 
this  is  not  always  possible.  Consider  a  D-Set  instance  generated  from  a  nonplansu'  graph  such 
that  each  edge  is  a  subset  containing  its  two  vertices  as  only  elements.  No  RB  instance  can 
be  mapped  to  such  a  D-Set  instance.  For  if  there  were  such  an  RB  instance,  the  geometric 
dual  of  the  planar  subdivision  it  defines  would  contain  a  planar  embedding  for  the  original 
nonplanar  graph.  This  is  an  impossibility,  hence  we  have  a  contradiction. 

Thus  the  set  of  RB  instances  constitutes  a  proper  subset  of  the  set  of  D-Set  instances;  in 
other  words,  RB  is  isomorphic  to  a  subproblem  of  D-Set.  Therefore,  the  NP-completeness 
of  RB  does  not  follow  directly  from  that  of  D-Set  established  earlier.  Fortunately,  however, 
D-Set  has  an  equivalent  subproblem  which  is  isomorphic  to  a  subproblem  of  RB  under 
That  isomorphism  provides  us  with  the  NP-completeness  of  Region  Basis. 

Theorem  3  Region  Basis  is  NP-complete. 

Proof.  That  RB  €  NP  is  easy  to  verify,  based  on  the  fact  mentioned  in  Section  3.1  that 
the  number  of  regions  in  the  planar  subdivision  is  at  most  quadratic  in  the  total  size  of  the 
polygons. 

Let  (X,  C)  be  a  D-Set  instance,  where 

X  =  {a:i,X2,...,a^m}; 

C  =  {Si,S2,‘“iSn}  Q  2^. 

Without  loss  of  generality,  we  make  two  assumptions 

\JSi  =  X  and  n  =  0, 

t=l  i=l 

because  elements  contained  in  none  of  the  subsets  or  contained  in  all  subsets  can  always 
be  removed  from  any  discriminating  set  of  {X,  C).  Now  add  in  a  new  element  a  ^  X  and 
consider  the  D-Set  instance  (A"  U  {a},C'),  where  C  =  {/S',  U  {a}  |  1  <  f  <  n}.  Clearly 
{X  U  {a},C7')  and  {X,C)  have  the  same  set  of  irreducible  discriminating  sets'*  and  hence 
they  are  considered  equivalent. 

The  planar  subdivision  defined  by  the  constructed  RB  instance  for  (A  U  {a},C")  taJces 
the  configuration  shown  in  Figure  9(a):  A  rectangular  region  is  divided  by  a  horizontal  line 
segment  into  two  identical  regions  of  which  the  bottom  region  is  named  u;(a);  the  top  region 
is  further  divided,  this  time  by  vertical  line  segments,  into  2m  —  1  identical  regions  of  which 
the  odd  numbered  ones,  from  left  to  right,  are  named  u;(xi), . . .  ,u?(xm)  respectively.  Remove 
those  m— 1  unnamed  regions  on  the  top.  For  1  <  i  <  n  define  polygon  Pi  to  be  the  boundary 
of  the  union  of  all  regions  w(x),  x  €  5,-  U  {a}.  It  should  be  clear  that  Pi  is  indeed  a  polygon; 
and  the  two  assumptions  guarantee  that  Fi, . . . , form  the  desired  subdivision.  Note  the 
subdivision  consists  of  m  -|- 1  rectangular  regions  and  4m  +  2  vertices.  All  can  be  computed 

discriminating  set  D  is  said  to  be  irreducible  if  no  subset  D'  C  D  can  be  a  discriminating  set. 
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Figure  9:  Two  reductions  from  Discriminating  Set  to  Region  Basis. 

in  time  0(m),  given  the  coordinates  of  the  four  vertices  of  the  bounding  rectangle.  Thus  the 
reduction  takes  time  @(mLi  |5,|). 

It  is  clear  that  C  has  a  discriminating  set  of  size  I  or  less  if  and  only  if  there  is  a  region 
basis  of  the  same  size  for  Pi, . . . ,  P„.  Hence  we  have  proved  the  NP-completeness  of  RB.  □ 

The  above  proof  implies  that  we  may  regard  Discriminating  Set  and  Region  Basis  as 
equivalent  problems.  Note  that  the  polygons  Pi,...,P„  in  Figure  9  are  not  convex;  will 
Region  Basis  become  P  when  all  the  polygons  are  convex?  This  question  is  answered  by  the 
following  corollary. 

Corollary  2  Region  Basis  remains  NP-complete  even  if  all  the  polygons  are  convex. 

Proof.  Same  as  the  proof  of  Theorem  3  except  that  we  use  the  planar  subdivision  shown 
in  Figure  9(b).  (The  vertices  of  the  subdivision  partition  an  imaginary  circle  (dotted)  into 
2n  equal  arcs.)  □ 

3.4  Approximation 

Sometimes  we  can  derive  a  polyi  amial-time  approximation  algorithm  for  the  NP-complete 
problem  at  hand  from  some  existing  approximation  algorithm  for  another  NP-complete  prob¬ 
lem,  by  reducing  one  problem  to  the  other.  In  fewer  cases,  where  the  reduction  preserves 
the  solutions,  namely,  every  instance  of  the  original  problem  and  its  reduced  instance  have 
the  same  set  of  solutions,  any  approximation  algorithm  for  the  reduced  problem  together 
with  the  reduction  will  solve  the  original  problem.  The  problem  to  which  we  will  reduce 
Discriminating  Set  is  Hitting  Set: 

Hitting  Set 

Given  a  collection  C  of  subsets  of  a  finite  set  X,  find  a  minimum  hitting  set  for  C,  i.e.,  a 
subset  H  C  X  of  minimum  cardinality  such  that  /f  n  5  0  for  all  5  €  C? 

Karp  [36]  shows  Hitting  Set  to  be  NP-complete  by  a  reduction  from  Vertex  Cover.  The 
reducibility  from  D-Set  to  Hitting  Set  follows  a  key  fact  we  observed  when  proving  Theorem  2: 
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The  intersections  of  a  finite  set  D  with  two  finite  sets  S  and  T  are  not  equal  if  and  only  if 
D  intersects  their  symmetric  difference  SAT.  Given  a  D-Set  instance,  the  corresponding 
Hitting  Set  instance  is  constructed  simply  by  replacing  all  the  subsets  with  their  pairwise 
symmetric  differences.  Thus  every  discriminating  set  of  the  original  D-Set  instance  is  also  a 
hitting  set  of  the  constructed  instance,  and  vice  versa. 

The  approximability  of  Hitting  Set  can  be  studied  through  another  problem,  Set  Covering: 

Set  Covering 

Given  a  collection  C  of  subsets  of  a  finite  set  X,  find  a  minimum  cover  for  X,  i.e.,  a  subcol¬ 
lection  C’  QC  of  minimum  size  such  that  Usec'  S  =  X? 

This  problem  is  also  shown  to  be  NP-complete  by  Karp  by  a  reduction  from  Exact  Cover 
by  3-Sets  [36].  A  greedy  approximation  algorithm  for  this  problem  due  to  Johnson  [34]  and 
Lovasz  [43]  guarantees  to  find  a  cover  C  for  X  with  ratio 

or  simply  <  In  |X| -|- 1, 

where  C*  is  a  minimum  cover  and  H(k)  =  Hk  =  7,  known  as  the  Arth  harmonic 

number.  The  algorithm  works  by  selecting,  at  each  stage,  a  subset  from  C  that  covers  the 
most  remaining  uncovered  elements  of  X.  We  refer  the  reader  to  [15]  for  a  general  analysis 
of  the  greedy  heuristic  for  Set  Covering. 

Hitting  Set  and  Set  Covering  are  duals  to  each  other — the  roles  of  set  and  element  in  one 
problem  just  get  switched  in  the  other.  More  specifically,  let  a  Hitting  Set  instance  consist 
of  some  finite  set  X  and  a  collection  C  of  its  subsets;  its  dual  Set  Covering  instaince  then 
consists  of  a  set  C  and  a  collection  of  its  subsets  X  where 

C  =  {5|5€C'}  and  X  =  {x\xeX}, 

and  where  each  subset  x  is  defined  as®’® 

X  =  {  5  1  S  €  C  and  5  9  x  }. 

Intuitively  speaJdng,  the  element  x  €  X  “hits”  the  subset  5"  €  C  in  the  original  instance  if 
and  only  if  the  subset  x  “covers”  the  element  S  €  C  in  the  dual  instance.  Thus  it  follows 
that  if  C  X  is  a  hitting  set  for  C  if  and  only  ifH={x\x€H}isa,  cover  for  C.  Hence  the 
corresponding  greedy  algorithm  for  Hitting  Set  selects  at  each  stage  an  element  that  “hits” 
the  most  remaining  subsets.  It  is  clear  that  the  approximation  ratio  for  Hitting  Set  becomes 
H{maXxex  |  {  S’  |  5  €  C  and  S  9  x  }  |)  or  In  jCj  -I- 1.^ 

'^According  to  this  definition,  x  =  y  may  hold  for  two  different  elements  x  ^  y.  In  this  case  only  one 
subset  is  included  in  X. 

®This  definition  also  establishes  the  duality  between  D-Set  and  a  known  NP-complete  problem  called 
Minimum  Test  Set  (see  [25]).  Given  a  collection  of  subsets  of  a  finite  set,  Minimum  Test  Set  asks  for  a 
mininum  subcollection  such  that  exactly  one  from  each  pair  of  distinct  elements  is  contained  in  some  subset 
from  this  subcollection. 

^Kolaitis  and  Thakur  [39]  syntactically  define  a  class  of  NP-complete  problems  with  logarithmic  approx¬ 
imation  algorithms  that  contains  Set  Covering  and  Hitting  Set,  and  show  that  Set  Covering  is  complete  for 
the  class. 
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As  a  shcurt  summary,  the  greedy  heuristic  on  a  Discriminating  Set  instance  {X^C)  works 
by  finding  a  hitting  set  for  the  instance  (X,  {S  AT  \  S,T  €C}).  Since  an  element  can  ap- 
pear  in  at  most  such  pairwise  symmetric  differences,  where  n  =  |C|,  the  approximation 
ratio  attained  by  this  heuristic  is  ln[^J  +  1  <  21nn.  The  same  ratio  is  attained  for  Region 
Basis  by  the  heuristic  that  selects  at  each  step  a  region  discriminating  the  most  remaining 
pairs  of  polygons,  where  n  is  now  the  number  of  polygons. 

The  greedy  algorithm  for  Set  Covering  (dually  for  Hitting  Set)  can  be  carefully  imple¬ 
mented  to  run  in  time  0(I2s€C  I'S'I)  [17].  The  reduction  from  a  D-Set  instance  (X,  C)  to  a 
Hitting  Set  instance  takes  time  0(|Cpmaxs£c  |5|).  Combining  the  time  complexity  of  the 
gemnetric  preprocessing  in  Section  3.1,  we  can  easily  verify  that  Region  Basis  can  be  solved 
in  time  0{nm^  -I-  n*m*)  =  0(n*m^),  where  n  and  m  are  the  number  and  size  of  polygons 
respectively. 

In  the  remainder  of  this  subsection  we  establish  the  hardness  of  approximating  D-Set  and 
hence  Region  Basis.  Both  problems  allow  the  same  approximation  ratio  since  the  reductions 
firom  one  to  another  do  not  change  the  number  of  subsets  (or  polygons)  in  an  instance.  First 
we  should  note  that  the  ratio  bound  H(maxscc  1^1)  of  the  greedy  algorithm  for  Set  Covering 
is  actually  tight;  an  example  that  makes  the  algorithm  achieve  this  ratio  for  arbitrarily  large 
maxsec  |<S'|  is  given  in  [34]. 

Next  we  present  a  reverse  reduction  from  Hitting  Set  to  D-Set  to  show  that  an  algorithm 
for  D-Set  with  approximation  ratio  clogn  can  be  used  to  obtain  an  algorithm  for  Hitting 
Set  with  ratio  clogn  -|-  O(loglogn),  where  c  >  0  is  any  constant  and  n  is  the  number  of 
subsets  in  an  instance.  Afterwards,  we  will  apply  some  recent  results  on  the  hardness  of 
approximating  Set  Covering  (and  thus  Hitting  Set). 

Lemma  2  For  any  c>  0,  if  c\ogn  is  the  approximation  ratio  of  Discriminating  Set,  then 
Hitting  Set  can  be  approximated  with  ratio  clogn  -)-  O(loglogn). 

Proof.  Suppose  there  exists  an  algorithm  A  for  D-Set  with  approximation  ratio  c  log  n.  Let 

(X,  C)  be  an  arbitrary  instance  of  Hitting  Set,  where  C  =  {Si, . . . ,  Sn}  C  2^ ,  and  let  n  =  jCj. 
To  construct  a  D-Set  instance,  we  first  make  /(n)  isomorphic  copies  (Xi,  Ci), . . . ,  (X/(„),  C/(„)) 
of  (X,^)  such  that  X,-  0  Xj  =  0  for  1  <  i  j  <  /(«)•  Here  /  is  an  as  yet  undetermined 
function  of  n  upper  bounded  by  some  polynomial  in  n.  Now  consider  the  enlarged  Hitting  Set 
instance  (X',C0  =  Every  hitting  set  W  of  (X',C")  has  W  =  u£f 

where  Hi  is  a  hitting  set  of  (X,-,  C,),  1  <  i  <  n;  so  from  H'  we  can  obtain  a  hitting  set  H  of 
(X, C)  with  |/f|  <  \H'\I f{n)  merely  by  taking  the  smallest  one  of  Hi,..., Hf(n). 

Next  we  introduce  a  set  A  consisting  of  new  elements  oi,  02, . . . ,  aiog(n/(n))  ^  X';  and  for 
1  <  i  <  n/(n)  define  auxiliary  sets  A,: 

A,  =  {  flj  1 1  <  y  <  log(n/(n))  and  the  jth  bit  of  the  binary  representation  of  t  —  1  is  1 }. 

It  is  not  hard  to  see  {ai,...,uiog(T./(n))}  must  be  a  subset  of  any  discriminating  set  for 
Ai,..., An/(n);  therefore  it  is  the  minimum  one  for  these  auxiliary  sets.  The  constructed 
D-Set  instance  is  then  defined  to  be  {X",C"),  where 

X"  =  Xi  U  UX/(„)  U  {01,02,..., aiog(„/(n))}; 

C"  =  {  T  U  A(, •  I  T  €  Ci  and  T  =  Sj}  U  {  Ai, . . . ,  A„/(„) }. 
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It  is  easy  to  verify  that  every  discriminating  set  of  (X",  C")  is  the  union  of  A  and  a  hitting 
setof  (X',C"). 

Now  run  algorithm  A  on  the  instance  {X'\  C")  and  let  D  be  the  discriminating  set  found. 
Then 

<  clog(|C"l)  =  clog(2n/(n)), 


where  D*  is  a  minimum  discriminating  set.  From  the  construction  of  (X",  C")  we  know 
that  D  =  H\\J  '"KJ  ■ff/(n)  U  A  and  D*  =  /fj  U  •  •  •  U  ^/(„)  U  A,  where  for  1  <  t  <  n,  /fi 
and  H*  are  some  hitting  set  and  some  minimum  hitting  set  of  (X,,C,),  respectively.  Let  Hk 
satisfy  \Hk\  =  min/^i^  and  thus  let  H  =  Hk  he  a.  hitting  set  of  (X,C7).  Also,  let  H*  with 

=  •  •  •  =  be  a  minimum  hitting  set  of  (X,C).  Then 


S/l.’  Itfil  +  \A\ 
iffri  + 1-41 

/(n)  •  1^1  +  log(n/(n)) 
/(n)-|^*|  +  log(n/(n))‘ 


Combining  the  two  inequalities  above  generates: 


m 


<  c\og{2nf{n))  + 


(clog(2n/(n))  -  1) » log(n/(n)) 


<  clog(2n/(n))  + 


fin) .  |tf-| 

(clog(2n/(n))  -  1)  •  log(n/(n)) 

/(«) 

=  clogn  +  [c  +  clog/(n)  +  +  ‘°8”  +  l°g/("))). 

I  /(n)  J 


Setting  /(n)  =  log^n,  all  terms  in  the  brackets  can  be  absorbed  into  0(log  log  n)  after  simple 
manipulations  on  asymptotics  [29];  thus  we  have 

<  c log  n  +  0(log  log  n). 

□ 


Though  Set  Covering  has  been  extensively  studied  since  the  mid  70’s,  essentially  nothing 
on  the  hardness  of  approximation  was  known  until  very  recently.  The  results  of  [4]  imply 
that  no  polynomial  approximation  scheme  exists  unless  P  =  NP.  Based  on  recent  results 
from  interactive  proof  systems  and  probabilistically  checkable  proofs  and  their  connection 
to  approximation,  several  asymptotic  improvements  on  the  hardness  of  approximating  Set 
Covering  have  been  made.  In  particular,  Lund  and  Yannakakis  [44]  showed  that  Set  Covering 
cannot  be  approximated  with  ratio  clogn  for  any  c  <  ^  unless  NP  C  DTIME(nP®*’'‘°*’‘);  Bel- 
lare  et  al.  [7]  showed  that  approximating  Set  Covering  within  any  constant  is  NP-complete, 
and  approximating  it  within  clogn  for  any  c  <  |  implies  NP  C  DTIME(n*°**°*’‘).  Based  on 
their  results  and  by  Lemma  2,  we  conclude  on  the  same  hardness  of  approximating  D-Set 
and  Region  Basis: 
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Theorem  4  Discriminating  Set  and  Region  Basis  cannot  be  approximated  by  a  polynomial¬ 
time  algorithm  with  ratio  bound  clog  n  /or  any  c  <  |  unless  NP  C  or /or 

anye<l  unless  NP  C  DTIME(n‘»«*-«’‘). 

Following  the  above  theorem,  the  ratio  21nn  »  1.39  log  n  of  the  greedy  algorithm  for 
D-Set  remains  asymptotically  optimal  if  NP  is  not  contained  in  DTIME(n**®'^*®*’‘). 


3.5  More  on  Discriminating  Set 

Now  let’s  come  back  to  where  we  left  the  discussion  on  the  subproblems  of  D-Set  in  Sec¬ 
tion  3.2;  it  has  not  been  settled  whether  D-Set  remains  NP-complete  when  every  subset  S 
in  the  collection  C  satisfies  |5|  <  2.  We  now  prove  that  this  subproblem  is  NP-complete. 

Here  we  look  at  a  special  case  of  this  subproblem,  namely,  a  "subsubproblem”  of  D-Set, 
subject  it  to  two  restrictions:  (1)  0  €  C7  and  (2)  |5|  =  2  for  all  nonempty  subsets  S  €  C. 
Let’s  call  this  special  case  0-8  D-Set.  If  0-2  D-Set  is  proven  to  be  NP-complete,  so  will  be 
the  original  subproblem. 

It’s  quite  intuitive  to  imderstand  a  0-2  D-Set  instance  in  terms  of  a  graph  G  =  (V,E), 
where  V  =  X,  the  finite  set  of  which  every  5  €  C  is  a  subset,  and 

E  =  {(u,t;)  I  {«,»}  e  C}. 

In  other  words,  each  element  of  the  set  X  corresponds  to  a  vertex  in  G  while  each  subset, 
except  0,  corresponds  to  an  edge.  Clearly  this  correspondence  from  all  0-2  D-Set  instances 
to  all  graphs  is  one-to-one.  Since  any  discriminating  set  D  for  C  has 

D  n  5  ^  D  n  0  =  0,  for  all  5  €  C  and  5/0, 

D  must  be  a  vertex  cover  for  G.  Let  d{u,v)  be  the  distance,  i.e.,  the  length  of  the  shortest 
path,  between  vertices  u,v  in  G  (or  oo  if  u  and  v  are  disconnected).  A  3-independent  set 
in  G  is  a  subset  I  QV  such  that  d(u,  v)  >  3  for  every  pair  u,  v  €  /.  The  following  lemma 
captures  the  dual  relationship  between  a  discriminating  set  for  C  and  a  3-independent  set 
in  G. 

Lemma  3  Let  X  be  a  finite  set  and  C  a  collection  of  0  and  two-element  subsets  of  X. 
Let  G  —  (X,  E)  be  a  graph  with  E  =  { (u,  v)  |  {u,  v}  G  G  } .  Then  a  subset  D  C  X  is  a 
discriminating  set  for  C  if  any  only  if  X\D  is  a  3-independent  set  in  G. 

Proof.  Let  D  be  a  discriminating  set  for  G.  Assume  there  exist  two  distinct  elements 
(vertices)  u,v  E  X\D  such  that  d(u,  u)  <  3.  We  immediately  have  (u,  v)  ^  E  since  D  must 
be  a  vertex  cover  in  G;  so  d(u,  v)  =  2.  Hence  there  is  a  third  vertex,  say  w,  that  is  connected 
to  both  u  and  v;  furthermore,  w  €  D  holds  since  the  edges  (u,  w)  and  (u,  w)  must  be  covered 
by  D.  But  now  we  have  jD  D  {u,u;}  =  D  H  {u,  w}  =  {u>},  a  contradiction  to  the  fact  that  D 
is  a  discriminating  set. 

Conversely,  suppose  X\D  is  a  3-independent  set  in  G,  for  some  subset  D  C  X.  Then 
D  must  be  a  vertex  cover.  Suppose  it  is  not  a  discriminating  set  for  G.  Then  there  exist 
two  distinct  subsets  5i,  52  €  G  such  that  DOSi  =  Df)  82  =  {tn},  for  some  to  G  5.  Writing 
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5*1  =  and  Sj  =  we  have  d(u,v)  =  2;  but  in  the  meantime  u,v  E  X  \  D.  A 

contradiction  again.  □ 

This  lemma  tells  us  that  the  NP-completeness  of  0-2  D-Set,  and  therefore  of  our  remaining 
open  subproblem  of  D-Set,  follows  if  we  can  show  the  NP-completeness  of  3- Independent  Set. 
3- Independent  Set  is  among  a  family  of  problems  defined,  for  all  integers  ib  >  0,  as  follows: 

ib-iNDEPENDENT  SET  (A:-IS) 

Given  a  graph  G  =  (V^E)  and  an  integer  0  <  /  <  1V|,  is  there  a  k-independent  set  of  size 
at  least  I,  that  is,  is  there  a  subset  I  QV  with  |/|  >  /  such  that  d{u,v)  >  k  for  every  pair 
UjV  £  II 

Thus  2-IS  is  the  familiar  NP-complete  Independent  Set  problem.  We  will  see  in  Appendix  B 
that  every  problem  in  this  family  for  which  A:  >  3  is  also  NP-complete.  To  avoid  too  much 
divergence  from  0-2  D-Set,  let’s  focus  on  3-IS  only  here. 

Lemma  4  S-Independent  Set  is  NP-complete. 

Proof.  It  is  trivial  that  3-IS  €  NP.  To  show  NP-hardness,  we  reduce  Independent  Set 
(2-IS)  to  3-IS.  Let  G  =  iV,E)  and  0  <  /  <  |V|  form  an  instance  of  Independent  Set.  A 
graph  G'  is  constructed  from  G  in  two  steps.  In  the  first  step,  we  introduce  a  “midvertex” 

otu,v  for  each  edge  (tt,v)  €  E,  and  replace  this  edge  with  two  edges  (u,a„,„)  and  In 

the  second  step,  an  edge  is  added  between  every  two  midvertices  that  are  adjacent  to  the 
same  original  vertex.  More  formally,  we  have  defined  G'  =  (V',  E')  where 

y  =  Vu{o«.„  I  («,«)€£}; 

E'  =  { (a«,„, «)  I  (u,  v)  €  }  U  { (a„,„,  |  (u,  u)  (u,  u;)  6  }. 

Two  observations  are  made  about  this  construction.  First,  it  has  the  property  that  d'(u,  u)  = 
d(u,  i;)-H  holds  for  any  pair  of  vertices  u,  t;  6  V,  where  d  and  d!  are  the  two  distance  functions 
in  G  and  G'  respectively.  This  equality  can  be  verified  by  contr2Mliction.  Next,  if  (u,v)  6  JS, 
then  any  two  midvertices  and  av,y  have 

d  (0!u^,Qfv,v)  ^  d  "b  d  {,OCu,vi^v,y)  —  2? 

d  u)  —  Ott.v)  "b  d^(Qu,v>  n)  ^  2j 

d  (Ov.yj  n)  =  d  (ofv,y,  0(u,v)  d  (Qy,y,u)  ^  2. 

Note  strict  ‘<’s  appear  in  above  three  inequalities  when  x  =  v  or  y  =  u,  and  in  the  first 
inequality  when  x  =  y.  It  is  not  difficult  to  see  that  the  whole  reduction  can  be  done  in  time 
0(|Vp).  Figure  10  illustrates  an  example  of  the  reduction. 

We  claim  that  G  has  an  independent  set  I  of  size  at  least  I  if  and  only  if  G'  has  a  3- 
independent  set  P  of  the  same  size.  Suppose  I  with  |/|  >  /  is  an  independent  set  in  G. 
Then  I  is  also  a  3-independent  set  in  G'.  This  follows  from  our  first  observation.  Conversely, 
suppose  r  with  |/'|  >  /  is  a  3-independent  set  in  G'.  Then  the  set  I,  produced  by  replacing 
eadi  midvertex  a„,„  G  /'  with  either  u  or  v,  is  an  independent  set  in  G.  To  see  this,  assume 
there  exists  two  vertices  u,t;  €  /  such  that  d(u,v)  =  1.  Thus  d'{u,v)  =  d{u,v)  +  1  =  2  <  3; 
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(•)  (b) 

Figure  10:  An  example  of  reduction  from  Independent  Set  to  3-IS.  (a)  An  instance  of  Independent 
Set.  (b)  The  constructed  3-IS  instance:  a  vertices  are  added  to  increase  the  distances  between  the 
original  vertices  by  exactly  one. 

so  either  u  or  v,  or  both,  must  have  replaced  some  midvertices  in  Let  s,  t  €  /'  be  the  two 
vertices  corresponding  to  u  and  v  before  the  replacement  respectively;  that  is,  s  =  u  or 
and  t  =  V  or  for  some  x,  y  €  V.  According  to  our  second  observation,  we  always  have 
<  2.  Thus  we  have  reached  a  contradiction,  since  s,f  6  That  |/'|  =  |/|  >  /  is  easy 
to  verify  in  a  similar  way.  □ 

Combining  Lemmas  3  and  4,  we  have  the  NP-completeness  of  0-2  D-Set;  this  immediately 
resolves  the  complexity  of  our  remaining  subproblem  of  D-Set: 

Theorem  5  D-Set  remains  NP-complete  even  */|5|  <  2  for  all  S  €  C. 


3.6  Experiments 

For  geometric  preprocessing,  we  implemented  the  plane  sweep  algorithm  by  Nievergelt  and 
Preparata  [47].  We  modified  the  original  algorithm  so  that  the  containing  polygons  of  each 
swept  region  are  maintained  and  propagated  along  during  the  sweeping.®  The  greedy  approx¬ 
imation  algorithm  for  Set  Covering  was  implemented  with  a  linked  list  to  attain  the  running 
time  0(Es€C  I-S'I)-  All  code  was  written  in  Common  Lisp  and  was  run  on  a  Sparcstation 
IPX. 

We  discuss  simulation  results  on  random  polygons.  These  simulations  empirically  study 
how  the  number  of  sampling  points  varies  with  the  “density”  of  polygons  in  the  plane. 
The  results  suggest  that  the  point  sampling  approach  is  most  effective  at  sensing  polygonal 
objects  that  have  highly  overlapping  poses.  Experiments  on  a  Zebra  robot  are  underway  and 
the  results  will  be  presented  in  the  near  future. 

*This  implementation  has  the  same  worst-case  running  time  as  a  different  version  described  in  Section  3.1 
which  obtains  the  containment  information  by  traversing  the  planar  subdivision  after  the  sweeping.  But  the 
implementation  version  is  usually  more  eflScient  in  practice. 
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3.6.1  Simulation  Results 


To  generate  random  polygons,  we  precomputed  an  arrangement  of  a  large  number  (such  as 
100)  of  random  lines  using  Edelsbrunner  and  Guibas’s  topological  sweeping  algorithm  [20]. 
A  random  polygon  was  extracted  as  the  first  “valid”  cycle  occurring  in  a  random  walk 
on  this  line  arrangement,  after  being  seeded  to  some  random  perimeter.  By  “valid”  we 
mean  that  the  number  of  vertices  in  the  cycle  was  no  less  than  some  small  random  integer. 
This  constraint  was  introduced  merely  to  allow  a  proper  distribution  of  polygons  of  various 
sizes,  for  otherwise  triangles  and  quadrilaterals  would  be  generated  with  high  probabilities 
according  to  our  observations.  In  a  sample  run,  a  group  of  1000  polygons  generated  (by  this 
method)  from  an  arrangement  of  100  random  lines  had  sizes  in  the  range  3-30,  with  mean 
5.545  and  standard  deviation  3.225. 

All  random  polygons  (or  all  random  poses  of  a  single  polygon)  in  a  test  were  bounded 
by  some  square,  so  that  the  “density”,  i.e.,  the  degree  of  overlap,  of  these  polygons  mainly 
depended  on  their  number  as  well  as  on  the  ratio  between  their  average  area  and  the  size  of  the 
boimding  square.  Since  polygons  were  generated  randomly,  the  average  area  could  be  viewed 
as  approximately  proportional  to  the  square  of  the  average  perimeter.  The  configuration  of 
each  polygon,  say  P,  was  assumed  to  obey  a  “uniform”  distribution  inside  the  square.  More 
specifically,  the  orientation  of  P  was  first  randomly  chosen  from  [0, 2ir);  the  position  of  P  was 
then  randomly  chosen  from  a  rectangle  inside  the  square  consisting  of  all  feasible  positions 
at  that  orientation.^ 

To  be  robust  against  sensor  noise,  the  sampling  point  of  every  region  in  the  region  basis 
was  selected  as  the  center  of  a  maximum  inscribed  circle  in  that  region.  In  other  words, 
this  sampling  point  had  the  maximum  distance  to  the  polygon  bounding  that  region.  It 
is  not  difficult  to  see  that  such  a  point  must  occur  at  a  vertex  of  the  generalized  Voronoi 
diagram  inside  the  polygon,  also  called  its  internal  skeleton  or  medial  axis  function}^  Also 
for  sensing  robustness,  regions  with  area  less  than  some  threshold  were  not  considered  at  the 
stage  of  region  basis  computation.^^  Though  this  thresholding  traded  off  the  completeness  of 
sampling,  it  almost  never  resulted  in  the  failure  of  finding  a  region  basis  once  the  threshold 
was  properly  set. 

The  first  two  groups  of  six  tests  gave  a  sense  of  the  number  of  sampling  points  required 
when  polygons  are  sparsely  distributed  in  the  plane.  The  results  are  summarized  in  Table  4. 
Every  test  in  group  (a)  was  conducted  on  distinct,  i.e.,  non-congruent,  random  polygons 
with  perimeters  between  |  and  |  of  the  width  of  the  bounding  square;  every  test  in  group 
(b)  was  conducted  on  distinct  poses  of  a  single  polygon  with  perimeter  equal  to  |  of  the 
width  of  the  square.  The  scene  of  the  la«t  test  from  each  group  is  displayed  in  Figure  11. 

Without  any  surprise,  the  number  of  sampling  points  found  were  around  half  of  the 

®If  the  diameter  of  P  is  greater  than  the  width  of  the  square,  then  not  every  orientation  is  necessarily 
feasible.  However,  this  situation  was  avoided  in  our  simulations. 

*°The  construction  of  the  internal  skeleton  of  a  polygon  is  a  special  case  of  the  construction  of  the  gen¬ 
eralized  Voronoi  diagram  for  a  set  of  line  segments,  for  which  0(n  log  n)  algorithms  are  given  in  [24],  [38], 
and  [56].  Since  the  maximum  region  size  for  a  region  basis  turned  out  to  be  very  small  in  the  simulations, 
we  only  implemented  an  O(n^)  brute  force  algorithm. 

^*We  thresholded  on  the  region  area  rather  than  the  radius  of  a  maximum  inscribed  circle  merely  to  avoid 
the  inefficient  computation  on  the  latter  for  all  the  regions  in  the  pluiar  subdivision. 
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#  polys 

#  regions 

#  sampling 
points 

50 

320 

31 

60 

500 

37 

70 

594 

41 

80 

783 

46 

90 

973 

47 

100 

1422 

51 

#  polys 

#  regions 

#  sampling 
points 

50 

362 

36 

60 

609 

34 

70 

741 

34 

80 

1061 

39 

90 

1125 

49 

100 

1643 

61 

(a)  (b) 


Table  4:  Tests  on  sampling  sparsely  distributed  random  polygons/poses.  The  twelve  tests  were 
divided  into  two  groups:  (a)  All  polygons  in  each  test  were  distinct,  with  perimeters  between  j 
and  I  times  the  width  of  the  bounding  square,  (b)  All  polygons  in  each  test  represented  distinct 
random  poses  of  a  same  polygon.  The  polygon  perimeter  was  uniformly  |  times  the  side  length  of 
the  square  for  all  six  tests  in  the  group. 


(a)  (b) 

Figure  11:  Sampling  100  sparsely  distributed  random  polygons/poses,  (a)  The  scene  of  the  last 
test  from  group  (a)  in  Table  4:  There  are  1422  regions  in  the  planar  subdivision  and  51  sampling 
points  (drawn  as  dots)  to  discriminate  the  100  polygons,  (b)  The  scene  of  the  last  test  from 
group  (b)  in  Table  4:  There  are  1643  regions  in  the  planar  subdivision  and  61  sampling  points  to 
discriminate  the  100  poses. 
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number  of  polygons,  for  all  twelve  tests  in  Table  4.  This  supports  the  fact  that,  for  n 
sparsely  distributed  polygons  in  the  plane,  the  minimum  number  of  sampling  points  turns 
out  to  be  0(n).  As  we  can  see  from  Figure  11,  in  such  a  situation  every  polygon  intersects 
at  most  a  few,  or  more  precisely,  a  constant  number  of,  other  polygons.  In  other  words,  the 
number  of  polygon  pairs  distinguishable  by  any  single  region  in  the  planar  subdivision  is 
0(.^);  but  there  axe  such  pairs  in  total!  Thus  sensing  by  point  sampling  is  inefficient  in 
a  situation  with  a  large  number  of  spau'sely  distributed  polygons. 

The  next  two  groups  of  six  tests  were  on  polygons  much  mor^  densely  distributed  in  the 
plane,  and  the  results  are  given  in  Table  5.  In  these  two  groups  of  tests,  we  used  a  bounding 


#  polys 

#  regions 

#  sampling 
points 

25 

776 

8 

30 

998 

8 

35 

1270 

12 

40 

1759 

12 

45 

2153 

il 

50 

2678 

13 

#  polys 

#  regions 

#  sampling 
points 

10 

264 

4 

20 

1121 

7 

25 

1781 

6 

30 

2796 

7 

35 

3655 

9 

40 

4995 

9 

(c)  (d) 

Table  5:  Tests  on  sampling  densely  distributed  random  polygons/poses,  divided  into  two  groups 
(c)  and  (d).  The  width  of  the  bounding  square  was  reduced  to  ^  times  the  width  of  the  square 
used  in  groups  (a)  and  (b)  in  Table  4.  In  group  (c)  all  polygons  in  each  test  were  distinct  with 
perimeters  in  the  range  between  ^  and  2  times  the  width  of  the  bounding  square.  In  group  (d)  all 
polygons  in  each  test  were  distinct  poses  of  the  same  polygon  ais  in  Figure  11(b). 

square  with  side  length  only  j  of  the  width  of  the  one  used  in  test  groups  (a)  and  (b).  Every 
test  in  group  (c)  was  conducted  on  distinct  polygons  with  perimeters  in  the  range  |-2  times 
the  side  length  of  the  bounding  square.  All  tests  in  group  (d)  were  distinct  poses  of  the  same 
polygon  used  in  the  last  test  of  group  (b).  Again  the  scene  of  the  last  test  from  each  group 
is  shown  in  Figure  12. 

All  twelve  tests  in  groups  (c)  <ind  (d)  except  the  last  one  in  group  (c)  found  sampling 
points  at  most  twice  the  lower  bound  flog  n] ,  while  the  first  test  in  group  (d)  found  exactly 
flogn]  sampling  points.  The  data  in  group  (d)  were  more  densely  distributed  than  the 
data  in  group  (c)  in  that  every  two  poses  intersected  each  other.  Since  an  extremely  dense 
distribution  of  polygons  may  cause  numerical  instabilities  in  the  plane  sweep  algorithm, 
smaller  numbers  of  polygons  were  tested  in  these  two  groups  than  were  tested  in  groups  (a) 
and  (b).  The  results  of  these  two  groups  of  tests  show  that  the  sampling  strategy  is  very 
applicable  to  sensing  densely  distributed  polygons. 
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Figure  12:  Sampling  densely  distributed  polygons.  The  bounding  square  has  width  i  times  the 
width  of  the  one  shown  in  Figure  11.  (c)  The  scene  of  the  last  test  from  group  (c)  in  Table  5: 
There  are  50  distinct  polygons  which  form  a  planar  subdivision  with  2678  regions,  and  which  can 
be  discriminated  by  13  sampling  points,  (d)  The  scene  of  the  last  test  from  group  (d)  in  Table  5: 
There  are  40  distinct  poses  of  the  polygon  from  Figure  11(b),  which  form  a  planar  subdivision  with 
4955  regions,  and  which  can  be  discriminated  by  9  sampling  points. 
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4  Conclusion 


The  analyses  and  experiments  in  this  paper  have  liiid  out  the  bases  for  two  general  sensing 
schemes  applicable  to  planar  objects  with  known  shapes.  The  first  scheme,  termed  sensing 
by  inscription,  determines  the  pose  of  an  object  by  finding  its  inscription  in  a  polygon 
of  geometric  constraints  derived  from  the  sensory  data.  An  implementation  may  use  a 
rotary  sensor  or  a  linear  CCD  array  combined  with  a  diverging  lens  to  obtain  the  necessary 
constraints.  In  particular,  two  supporting  cones  are  often  enough  to  detect  the  real  pose  of 
a  polygonal  object.  In  real  situations,  if  two  (or  more)  possible  poses  arise  from  a  two-cone 
inscription,  they  can  be  distinguished  by  point  sampling.^^ 

Though  only  the  inscription  of  a  convex  polygon  is  treated  the  extensions  to  any  cirbitrary 
polygon  and  any  polyhedron  should  be  straightforward.  In  the  second  case  the  3-D  cones 
are  defined  by  visible  vertices  as  well  as  edges.  But  the  extension  to  a  closed  and  piecewise 
smooth  curve  needs  further  study.  The  technique  Ccin  also  be  applied  in  object  recognition: 
A  finite  set  of  polygons  is  generally  distinguishable  by  inscription. 

Robustness  of  inscription  can  be  realized  by  allowing  some  tolerance  for  intersecting  locus 
curves.  That  is,  an  intersection  of  two  locus  curves  is  considered  to  be  a  real  pose  if  it  is 
attained  on  these  curves  at  orientations  which  differ  by  an  amount  less  than  the  tolerance. 

Future  work  will  involve  the  design  of  specialized  cone  sensors  or  other  sensors  suited  for 
inscription,  as  well  as  an  investigation  of  a  theoretical  framework  for  incorporating  sensing 
uncertainties  into  the  inscription  algorithms. 

The  second  scheme,  termed  sensing  by  point  sampling,  distinguishes  a  finite  set  of  poses 
by  examining  the  containment  of  a  number  of  points.  It  can  be  generalized  to  work  for 
objects  with  more  general  boundaries  as  long  as  the  planar  subdivision  can  be  efficiently 
constructed.  Robustness  mainly  depends  on  the  invariance  of  the  topology  of  the  planar 
subdivision.  Therefore  it  is  often  better  to  consider  only  regions  whose  maximum  inscribed 
circles  are  larger  than  a  certain  threshold. 

Future  work  will  combine  existing  mechanic2d  orienting  methods  with  point  light  detectors 
or  mechanical  probes  for  pose  detection. 
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(C) 


Figure  13:  A  worst  case  example  of  inscription.  The  three  half- planes  form  an  equilateral  triangle 
by  intersection  and  the  convex  n-gon  P  has  vertex  po  at  the  center  of  an  arc  equally  divided 
by  the  remaining  vertices  pi, . . .  iPn-i  of  P.  The  conditions  derived  in  this  appendix  on  radius  r 
and  measure  a  of  the  arc  guarantee  6n  different  poses  in  which  P  is  inscribed  in  (a) 

6(n  - 1)  poses  of  P  result  from  rotations  about  po  positioned  at  vertices  qo>  9i>  and  q2  respectively; 
and  (b)  6  other  poses  result  when  po,  pi  and  p„_i  are  incident  on  different  edges  of  Aqo9i92-  (c) 
illustrates  a  pose  of  the  first  type  such  that  vertex  p„  2  <  »  <  n  -  2,  is  on  edge  qiq2  during  the 
rotation  about  qo. 


Appendices 


A  A  Worst  Case  of  Inscription 

In  this  appendix  we  construct  an  example  where  a  convex  n-gon  embedded  in  three  half¬ 
planes  can  attain  6n  possible  poses. 

Symmetry  turns  out  to  play  a  central  role  in  the  construction.  Let  the  three  half-planes 
form  an  equilateral  triangle  Aqo9i93  by  intersection;  and  let  F  be  a  convex  polygon  with 
vertices  pq,pi, . . .  ,p„_i  in  counterclockwise  order  such  that  po  is  the  center  of  an  circ  of 
radius  r  and  measure  a,  and  pi, . . .  ,Pn-i  together  divide  this  arc  into  n  —  2  equal  pieces. 
The  reason  for  choosing  this  particular  shape  of  P  is  that  we  expect  to  obtain  2(n  —  1)  poses 
by  positioning  po  at  each  of  qo,  qi,  and  qj  and  rotating  P  about  po  inside  Aqoqiq2  such  that 
Pi, . . .  ,Pn-i  will  all  touch  the  opposite  edge  of  Aqo9i92  ex2w:tly  twice  during  the  rotation. 
Figure  13(a)  illustrates  the  first  of  the  sequence  of  2(n  —  1)  poses  resulting  from  the  rotation 
around  qo*  By  symmetry  we  already  have  6(n  —  1)  poses  in  total,  called  poses  of  the  first 
type.  The  remaining  6  poses  of  the  seconu  type  are  symmetric  to  each  other,  attained  when 
vertices  po,  pi  and  p„_i  are  on  different  edges  of  Aqoqiq2.  (See  Figure  13(b).)  So  our  task 
is  to  derive  the  conditions  on  radius  r  and  measure  a  of  the  arc  that  realize  these  two  types 
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of  poses. 

Let  us  first  look  at  poses  of  the  first  type.  Let  cf  >  0  be  the  altitude  of  Aqoqiq2  and 
0  an  acute  angle  with  cosP  =  and  suppose  po  coincides  with  qo.  We  require  d  <  r  and 
a<a  +  )?<|so  that  p\  in  the  initial  pose  is  not  on  edge  qoqx  in  order  to  allow  poses  of 
the  second  type  shown  in  Figure  13(b),  and  p„_i  lies  to  the  left  of  the  midpoint  on  edge  qiq2 
in  order  to  be  incident  on  this  edge  the  second  time  as  P  rotates  counterclockwise.  We  thus 
have  the  following  conditions  on  a  and  r: 


0<a<- 


cos(—  —  a)  <  -  <  1. 
6  r 


With  the  above  constraints,  we  still  need  to  make  sure  that  P  is  indeed  above  edge  qiq2 
when  each  crtex  p,-,  1  <  t  <  n  —  1,  becomes  incident  on  this  edge  during  the  rotation.  In 
fact,  it  suffices  to  make  sure  that  vertices  p,_i  and  p^+i  jure  above  edge  qiq2.  Figure  13(c) 
illustrates  the  case  in  which  p,-  lies  to  the  left  of  the  midpoint  on  edge  q\q2.  It  is  clear  that 
Pi-i  is  above  edge  q\q2\  Pi+i  is  above  edge  qiq2  if  and  only  if  <  Zpi9oPi+i  =  and  thus 
if  and  only  if 

d  a 

r  2(n  -  2) 

Starting  from  the  initial  pose  shown  in  Figure  13(a),  slide  po  down  along  edge  9092  and 
Pn-i  left  along  edge  9192  until  pi  touches  edge  9091;  then  we  get  the  pose  in  Figure  13(b). 
Observe  that  po  is  on  edge  9092  but  closer  to  90  than  to  92,  thereby  ensuring  the  other  5 
poses  of  the  second  type  by  symmetry. 

Hence  P  has  fin  feasible  poses  when  inscribed  in  A909192  if  the  following  two  conditions 
are  satisfied: 

0  <  a  <  -; 

0 

,  .  /  d  ^  \ 

d<r  <mm( — — - r, - - — I. 

^cos(l-a) 

In  fact,  the  above  conditions  are  also  necessary  for  the  given  shapes  of  P  Jind  Aqoqiq2- 

Using  the  above  idea,  we  can  similarly  construct  a  worst  case  exaimple  for  the  special 
case  of  inscription  in  which  two  of  the  three  bounding  lines  are  parallel;  it  was  stated  in 
Section  2.2.3  that  there  exist  up  to  4n  poses  for  a  convex  n-gon.  Let  the  three  bounding 
lines  be  li,l2  and  I3  such  that  /i  ||  I2  and  /a  ±  hihi  f^t  a  convex  n-gon  P  take  the  same 
shape  just  described.  Similarly  we  have  that  P  has  4n  feasible  poses,  4(n  —  1)  of  which  result 
from  rotations  about  po  incident  on  l\  and  l2  respectively,  and  4  of  which  result  when  pi 
on  one  of  l\  and  l2  and  p„_i  on  the  other,  for  any  radius  r  and  measure  a  that  satisfy  the 
following  two  conditions: 


<  r  <  mm{ 


2(n-2) 


where  |  €  [0,  |)  is  the  root  of  equation  2sin  |  =  cos  ,  and  d  is  the  distance  between  lx 
and  I2. 
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B  A;-Independent  Sets 

We  extend  Lemma  4  to  all  k-lS  with  it  >  3:  They  are  NP-complete  as  well.  The  proof 
we  will  present  is  indeed  a  generalization  of  the  proof  of  Lemma  4;  it  will  again  construct 
a  ib-IS  instance  with  graph  G*  from  an  instance  of  Independent  Set  with  graph  G  by  local 
^  replacement.  In  the  proof,  each  vertex  u  in  G  will  be  replaced  by  a  simple  path  Vv  of  fixed 

length  (depending  only  on  k)  that  has  v  in  the  middle  and  an  equal  number  of  auxiliary 
vertices  on  each  side;  and  each  edge  (u,  u)  will  be  replaced  by  four  edges  connecting  the  two 
►  end  vertices  on  Vu  with  the  two  end  vertices  on  Vvt  either  directly  or  through  a  “midvertex”. 

More  intuitively  speaking,  all  shortest  paths  between  pairs  of  vertices  in  G,  if  they  exist,  get 
elongated  in  G'  to  such  a  degree  that  (1)  (u,v)  is  an  edge  in  G  if  and  only  if  the  distance 
between  vertices  u  and  v  in  G'  is  less  than  k;  (2)  any  two  vertices  u'  and  v'  in  G'  with  a 
distance  of  at  least  k  can  be  easily  mapped  to  two  nonadjacent  vertices  in  G.  The  first 
condition  ensures  that  any  given  independent  set  in  G  will  be  a  Ar-independent  set  in  G', 
while  the  second  condition  ensures  the  construction  of  aa  independent  set  in  G  from  any 
given  ^-independent  set  in  G'. 

Lemma  5  k-Independent  Set  is  NP-complete  for  all  integers  k  >  Z. 

Proof.  Given  an  instance  of  Independent  Set  as  a  graph  G  =  (V,  E)  and  a  positive  integer 
^  <  IV^Ij  a  A:-IS  instance  is  constructed  by  two  consecutive  substitutions.  A  path 

p  Vk-3  • .  •  V1VV2 . . .  Vk-2,  if  k  even; 

"  (  Vk.4 . . .  V1VV2 .  • .  Vk-3,  if  k  odd, 

first  substitutes  for  vertex  u  €  V,  where  Vi, . . .  ,^*-3  (and  Vk-2  when  k  is  even)  are  auxiliary 
vertices.  And  then  a  set  of  four  edges 

{(«ife-3>  Vk-3),  {Uk-3,  Vk-2), 

{uk-2,  Vk-3),  («jt-2,  wjfc-2)},  if  k  even; 

4j  (Wfc— 3j  Q^u,v)> 

(t;fe_4,  a„,„),  (vjfe_3,  Ou.w)},  if  k  odd, 

substitute  for  each  edge  («,  u)  €  E,  where  is  an  introduced  midvertex.  Figure  14  shows 
two  subgraphs  after  applying  the  above  substitutions  on  edge  (u,  v)  €  E,  for  k  even  and  odd 
respectively. 

*  We  can  easily  verify  that,  for  any  pair  of  vertices  x  on  and  y  on  Vv,  both  k  even  and 

odd,  we  have  d'{x,  y)  <  d\u,  w)  =  A:  —  1  <  A:  if  (u,  v)  €  jE,  where  d!  is  the  distance  function 
defined  on  G'.  On  the  other  hand,  if  (u,v)  ^  E,  we  have  d'{u,v)  >  k  when  k  is  even  and 
■*  cf{u,v)>k  +  l  when  k  is  odd.  Thus  an  independent  set  /  in  G  is  also  a  fc-independent  set  in 

G'.  Conversely,  suppose  P  with  |/'|  >  /  is  a  ^-independent  set  in  G.  We  substitute  u  G  V  for 
every  auxiliary  vertex  «,  €  /  on  path  Vu,  and  u  or  v  for  every  mid  vertex  6  I  when  k  is 

odd.  Let  I  be  the  set  after  this  substitution.  It  needs  to  be  shown  that  I  is  an  independent 
set  in  G  and  j/j  =  |/'|  >  /.  This  is  obvious  for  the  case  that  k  is  even.  When  k  is  odd. 
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(!)  (b) 

Figure  14:  Two  subgraphs  resulting  from  the  described  substitutions  performed  on  edge  (u,  v)  £  E 
for  the  cases  that  (a)  k  is  even;  and  (b)  k  is  odd. 

however,  the  situation  is  a  bit  more  complicated  due  to  the  possible  occurrences  of  those  a 
vertices  in  We  observe,  for  smy  a^.u,  Qfu'.u'j  and  x  on  path  where  u,v,u’,v',w  €  V, 

Jb  —  3 

<^{oiu,vjx)  <  — - - h  3  <  A:,  if  tu  =  ti  or  u,  or  (u,  w)  €  E,  or  (v,  m)  €  E; 

d'(Q!„,v» <  4  <  A:,  if  (u,u')  6  E. 

In  fact,  these  two  conditions  guarantee  that  /  is  an  independent  set  in  G  and  |/|  =  |/'|, 
which  we  leave  for  the  reader  to  verify. 

The  reduction  can  be  done  in  time  0{k\V\  +  |£^|),  which  reduces  to  0(1^1  +  l-El)  if 
k  is  treated  as  a  constant,  in  contrast  to  the  time  0(|Kp)  required  for  the  reduction  from 
Independent  Set  to  3-IS.  This  time  reduction  is  due  to  the  fact  that  midvertices  corresponding 
to  the  same  vertex  in  V  no  longer  have  edges  between  each  other.  □ 

Since  1-IS  can  be  easily  solved  by  comparing  |V|  and  /,  we  are  now  ready  to  sum  up  the 
complexity  results  on  this  family  of  problems  in  the  following  theorem. 

Theorem  6  k-Independent  Set  is  in  P  if  k  =  l  and  NP-complete  for  all  k  >2. 
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